Quantum weights of dyons and of instantons with non-trivial holonomy 
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We calculate exactly functional determinants for quantum oscillations about periodic instantons 
with non-trivial value of the Polyakov line at spatial infinity. Hence, we find the weight or the 
probability with which calorons with non-trivial holonomy occur in the Yang-Mills partition func- 
tion. The weight depends on the value of the holonomy, the temperature, Aqcd, and the separation 
between the BPS monopoles (or dyons) which constitute the periodic instanton. At large separa- 
tion between constituent dyons, the quantum measure factorizes into a product of individual dyon 
measures, times a definite interaction energy. We present an argument that at temperatures below 
a critical one related to Aqcd , trivial holonomy is unstable, and that calorons "ionize" into separate 
dyons. 

PACS numbers: 11.15.-q,11.10.Wx,11.15.Tk 



There are two known generalizations of the standard self-dual instantons to non-zero temperatures. One is the 
periodic instanton of Harrington and Shepard [1] studied in detail by Gross, Pisarski and YafFe [2]. These periodic 
instantons, also called calorons, are said to have trivial holonomy at spatial infinity. It means that the Polyakov line 



assumes values belonging to the group center Z{N) for the SU{N) gauge group [34]. The vacuum made of those 
instantons has been investigated, using the variational principle, in ref. [3]. 

The other generalization has been constructed a few years ago by Kraan and van Baal [4] and Lee and Lu [5] ; it has 
been named caloron with non-trivial holonomy as the Polyakov line for this configuration does not belong to the group 
center. We shall call it for short the KvBLL caloron. It is also a periodic self-dual solution of the Yang-Mills equations 
of motion with an integer topological charge. In the limiting case when the KvBLL caloron is characterized by trivial 
holonomy, it is reduced to the Harrington-Shepard caloron. The fascinating feature of the KvBLL construction is 
that a caloron with a unit topological charge can be viewed as "made of" Bogomolnyi-Prasad-Sommerfeld (BPS) 
monopoles or dyons [6, 7]. 

Dyons are self-dual solutions of the Yang-Mills equations of motion with static {i.e. time- independent) action 
density, which have both the magnetic and electric field at infinity decaying as Therefore these objects carry 

both electric and magnetic charges (prompting their name). In the 3-|-l-dimensional SU{2) gauge theory there are in 
fact two types of self-dual dyons [8]: M and L with (electric, magnetic) charges (-I-, +) and (— , — ), and two types of 
anti-self-dual dyons AI and L with charges (+,— ) and (—,+), respectively. Their explicit fields can be found e.g. in 
ref. [9]. In the SU{N) theory there are 2N different dyons [8, 10]: Afi, Af2, ...Mn-^i ones with charges counted with 
respect to A^ — 1 Cartan generators and one L dyon with charges compensating those of Mi...Mn^i to zero, and 
their anti-self-dual counterparts. 

Speaking of dyons one implies that the Euclidean space-time is compactified in the 'time' direction whose inverse 
circumference is temperature T, with the usual periodic boundary conditions for boson fields. However, the 
temperature may go to zero, in which case the Ad Euclidean invariance is restored. 

Dyons' essence is that the A4 component of the dyon field tends to a constant value at spatial infinity. This 
constant A4 can be eliminated by a time-dependent gauge transformation. However then the fields violate the 
periodic boundary conditions, unless A4 has quantized values corresponding to trivial holonomy, i.e. unless the 
Polyakov line belongs to the group center. Therefore, in a general case one implies that dyons have a non-zero value 



I. MOTIVATION AND THE MAIN RESULT 
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FIG. 1: The action density of the KvBLL caloron as function of z,t at fixed x — y — 0, with the asymptotic value of A4, at 
spatial infinity v — 0.9nT, v — I.IttT. It is periodic in t direction. At large dyon separation the density becomes static (left, 
ri2 = 1.5/T). As the separation decreases the action density becomes more like a 4d lump (right, ri2 = 0.6/T). In both plots 
the L,M dyons are centered at zl = —Yri2/2nT, zm ~ vri2/2nT, xl,m ~ yL,M = 0. The axes are in units of temperature T. 



of at spatial infinity and a non-trivial holonomy. 



The KvBLL caloron of the SU{2) gauge group (to which we restrict ourselves in this paper) with a unit topological 
charge is "made of one L and one M dyon, with total zero electric and magnetic charges. Although the action 
density of isolated L and M dyons does not depend on time, their combination in the KvBLL solution is generally 
non-static: the L,M "constituents" show up not as 3d but rather as Ad lumps, see Fig. 1. When the temperature 
goes to zero while the separation between dyons remain fixed, these lumps merge, and the KvBLL caloron is reduced 
to the usual Belavin-Polyakov-Schwarz-Tyutin instanton [11] (as is the standard Harrington-Shepard caloron), plus 
corrections of the order of T. However, the holonomy remains fixed and non-trivial at spatial infinity. 

There is a strong argument against the presence of either dyons or KvBLL calorons in the Yang-Mills partition 
function at nonzero temperatures [2]. The point is, the 1-loop effective action obtained from integrating out fast 
varying fields where one keeps all powers of A4 but expands in (covariant) derivatives of A4 has the form [12] 

S^-^°°P = J d^x [PiAi) + E^fEiAi) + B^feiAi) + higher derivative terms] , 



PiA4) = ^v^(2.T-v)^ 



V = V'Ipf [for the SU{2) group] (2) 

mod 2ttT 



where the perturbative potential energy term ^(^4) has been known for a long time [2, 14], see Fig. 2. As follows 
from eq.(l) the trace of the Polyakov line is related to v as 

^TrL^cos^. (3) 

The zeros of the potential energy correspond to ^TrL = ±1, i.e. to the trivial holonomy. If a dyon has v 7^ 2TiTn 
at spatial infinity the potential energy is positive-definite and proportional to the id volume. Therefore, dyons and 
KvBLL calorons with non-trivial holonomy seem to be strictly forbidden: quantum fluctuations about them have an 
unacceptably large action. 



Meanwhile, precisely these objects determine the physics of the supersymmetric YM theory where in addition to 
gluons there are gluinos, i.e. Majorana (or Weyl) fermions in the adjoint representation. Because of supersymmetry, 
the boson and fermion determinants about L, M dyons cancel exactly, so that the perturbative potential energy (2) is 
identically zero for all temperatures, actually in all loops. Therefore, in the supersymmetric theory dyons are openly 
allowed. [To be more precise, the cancellation occurs when periodic conditions for gluinos are imposed, so it is the 
compactification in one (time) direction that is implied, rather than physical temperature which requires antiperiodic 
fermions.] Moreover, it turns out [10] that dyons generate a non-perturbative potential having a minimum at v = ttT, 
i.e. where the perturbative potential would have the maximum. This value of A4 corresponds to the holonomy 
TrL = at spatial infinity, which is the "most non-trivial"; as a matter of fact <TrL>= is one of the confinement's 
requirements. In the supersymmetric YM theory configurations having TrL = at infinity are not only allowed but 
dynamically preferred as compared to those with ^Tri = ±1. In non-supersymmetric theory it looks as if it is the 
opposite. 
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Nevertheless, it has been argued in ref. [15] that the perturbative potential energy (2) which forbids individual 
dyons in the pure YM theory might be overruled by non-perturbative contributions of an ensemble of dyons. For fixed 
dyon density, their number is proportional to the 3d volume and hence the non-perturbative dyon-induced potential 
as function of the holonomy (or of A4 at spatial infinity) is also proportional to the volume. It may be that at 
temperatures below some critical one the non-perturbative potential wins over the perturbative one so that the system 
prefcas < TrL >= 0. This scenario could then serve as a microscopic mechanism of the confinement-deconfinement 
phase transition [15]. It should be noted that the KvBLL calorons and dyons seem to be observed in lattice 
simulations below the phase transition temperature [16-18]. 

To study this possible scenario quantitatively, one first needs to find out the quantum weight of dyons or the 
probability with which they appear in the Yang-Mills partition function. Unfortunately, the single-dyon measure is 
not well defined: it is too badly divergent in the infrared region owing to the weak (Coulomb-like) decrease of the 
fields. What makes sense and is finite, is the quantum determinant for small oscillations about the KvBLL caloron 
which possesses zero net electric and magnetic charges. To find this determinant is the primary objective of this 
study. The KvBLL measure depends on the asymptotic value of A4 (or on the holonomy through eq.(3)), on the 
temperature T, on A, the scale parameter obtained through the renormalization of the charge, and on the dyon 
separation ri2- At large separations between constituent L,M dyons of the caloron, one gets their weights and their 
interaction. 

The problem of computing the effect of quantum fluctuations about a caloron with non-trivial holonomy is of the 
same kind as that for ordinary instantons (solved by 't Hooft [19]) and for the standard Harrington-Shepard caloron 
(solved by Gross, Pisarski and Yaffe [2]) being, however, technically much more difficult. The zero-temperature 
instanton is 0(4) symmetric, and the Harrington-Shepard caloron is 0(3) symmetric, which helps. The KvBLL 
caloron has no such symmetry as obvious from Fig. 1. Nevertheless, we have managed to find the small-oscillation 
determinant exactly. It becomes possible because we are able to construct the exact propagator of spin-0, isospin-1 
field in the KvBLL background, which by itself is some achievement. 

As it is well known [2, 19], the calculation of the quantum weight of a Euclidean pseudoparticle consists of three 
steps: i) calculation of the metric of the moduli space or, in other words, computing the Jacobian composed of zero 
modes, needed to write down the pseudoparticle measure in terms of its collective coordinates, ii) calculation of the 
functional determinant for non-zero modes of small fluctuations about a pseudoparticle, iii) calculation of the ghost 
determinant resulting from background gauge fixing in the previous step. Problem i) has been actually solved already 
by Kraan and van Baal [4]. Problem ii) is reduced to iii) in the self-dual background field [20] since for such fields 
Det{Wf^u) = Det(— D^)'*, where W^iy is the quadratic form for spin-1, isospin-1 quantum fluctuations and is the 
covariant Laplace operator for spin-0, isospin-1 ghost fields. Symbolically, one can write 

KvBLL measure = J ^(collective coordinates) • Jacobian • Det~ ' (;^^^) . Det(— 2?^) (4) 

where the product of the last two factors is simply Det~^(— Z?^) in the self-dual background. As usually, the 
functional determinants are normalized to free ones (with zero background fields) and UV regularized (we use the 
standard Pauli-Villars method). Thus, to find the quantum weight of the KvBLL caloron only the ghost determinant 

has to be computed. 

To that end, we follow Zarembo [21] and find the derivative of this determinant with respect to the holonomy or, 
more precisely, to v = y^]4fAf 1 The derivative is expressed through the Green function of the ghost field in 
the caloron background. If a self-dual field is written in terms of the Atiah-Drinfeld-Hitchin-Manin construction, 
and in the KvBLL case it basically is [4, 5], the Green function is generally known [23]- [25] and we build it explicitly 
for the KvBLL case. Therefore, we are able to find the derivative 9Det(— £>^)/9v. Next, we reconstruct the full 
determinant by integrating over v using the determinant for the trivial holonomy [2] as a boundary condition. This 
determinant at v = is still a non-trivial function of ri2 and the fact that we match it from the v ^ side is a serious 
check. Actually we need only one overall constant factor from ref. [2] in order to restore the full determinant at 
V 7^ 0, and we make a minor improvement of the Gross-Pisarski- Yaffe calculation as we have computed the needed 
constant analytically. 

Although all the above steps can be performed explicitly, at some point the equations become extremely lengthy 
typical expressions are several Mbytes long and so far we have not managed to simplify them such that they would 
fit into a paper. However, we are able to obtain compact analytical expressions in the physically interesting case of 
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FIG. 2: Potential energy as function of v/T. Two minima correspond to ^TrL — ±1, the maximum corresponds to TrL = 0. 
The range of the holonomy where dyons experience repulsion is shown in dashing. 



large separation between dyons, ri2 ^ 1/T. We have also used the exact formulae to check numerically some of the 
intermediate formulae, in particular at v — s- 0. 

If the separation is large in the temperature scale, ri2 3> 1/T, the final result for the quantum measure of the KvBLL 
caloron can be written down in terms of the 3d positions of the two constituent L,M dyons £1^2, their separation 
''12 = \zi — Z2I, the asymptotic of A4 at spatial infinity denoted by v S [0, 27rT] and v = 2TrT — v £ [0, 27rT], see 
eq.(80). We give here a simpler expression obtained in the limit when the separation between dyons is much larger 
than their core sizes: 



X exp[-27rri2P"(v)] exp [-y(3)p(^)l ^ (5) 

where the overall factor C is a combination of universal constants; numerically C = 1.031419972084. A is the scale 
parameter in the Pauli-Villars scheme; the factor g^^ is not renormalized at the one-loop level. 



Since the caloron field has a constant A4 component at spatial infinity, it is suppressed by the same 
perturbative potential P(v) as given by eq.(2). Its second derivative with respect to v is P"{v) = 

v — ttT^I — v — ttt(i + ^^ . If V is in the ranges between and ttT ^1 — -^^ or between ttT ^1 + ^^ and 

2ttT (corresponding to the holonomy not too far from the trivial, 0.787597< jlTrLj < 1) the second derivative P"{v) 
is positive, and the L and M dyons experience a linear attractive potential. Integration over the separation ri2 of 
dyons inside a caloron converges. We perform this integration in section VII, estimate the free energy of the caloron 
gas and conclude that trivial holonomy (v — 0, 27rT) may be unstable, despite the perturbative potential energy P(v). 

In the complementary range ttT < v < ttT + (or ^iTrL] < 0.787597), P"(v) is negative (see Fig. 2), 

and dyons experience a strong linear-rising repulsion. It means that for these values of v, integration over the dyon 
separations diverges: calorons with holonomy far from trivial "ionize" into separate dyons. 



II. THE KVBLL CALORON SOLUTION 



Although the construction of the self-dual solution with non-trivial holonomy has been fully performed indepen- 
dently by Kraan and van Baal [4] and Lee and Lu [5] we have found it more convenient for our purposes to use the 
gauge convention and the formalism of Kraan and van Baal (KvB) whose notations we follow in this paper. 

The key quantities characterizing the KvBLL solution for a general SU{N) gauge group are the A^— 1 gauge-invariant 
eigenvalues of the Polyakov line (1) at spatial infinity. For the SU{2) gauge group to which we restrict ourselves in 
this paper, it is just one quantity, e.g. TrL, eq.(3). In a gauge where A4 is static and diagonal at spatial infinity, 
i.e. ^4|^_^oo ~ it t^is asymptotic value v which characterizes the caloron solution in the first place. We 

shall also use the complementary quantity v = 27rT — v. Their relation to parameters lu,lu introduced by KvB [4] is 
LU = ijj — 4^ — ^ — LU. Both V and v vary from to 27rT. At v = 0, 2ttT the holonomy is said to be 'trivial', 
and the KvBLL caloron reduces to that of Harrington and Shepard [1]. 
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There are, of course, many ways to parametrize the caloron solution. Keeping in mind that we shaU be mostly inter- 
ested in the case of widely separated dyon constituents, we shall parametrize the solution in terms of the coordinates 
of the dyons' 'centers' (we call constituent dyons L and M according to the classification in ref. [9]): 

L dyon : Zi = 



M dyon : Z2 



T 



T ' 

dyon separation : Z2 — Zi = ri2, |ri2| = ttT p^, 

where p is the parameter used by KvB; it becomes the size of the instanton at v — > 0. We introduce the distances 
from the 'observation point' x to the dyon centers, 

f=x — zi=x + 2wri2, r = \r\, 

s = x — Z2=x — 2u)ri2, s = \s\. (6) 

Henceforth we measure all dimensional quantities in units of temperature for brevity and restore T explicitly only in 
the final results. 

The KvBLL caloron field in the fundamental representation is [4] (we choose the separation between dyons to be 
in the third spatial direction, ri2 = ''1263): 

= 5^,4 IV Y + ^€.^^30. ln$ + Re [(77I, - zr7^J(ri + iT2){d, + iy6,,4)x] , (7) 

where are Pauli matrices, r^J^^ are 't Hooft's symbols [19] with fjfj = Caij and fjl^, = —f]f;4 = iSav "Re" means 
2Re(M^) = W + and the functions used are 

tj) = — cos(27ra;4) + chch H shsh , ip = tp -\ — i^shsh -\ shch H shch , 

2rs rs s r 

ri2 f^-2^^^^si^ , sh\ ^ _ V 



^ V s r J ^ 

We have introduced short-hand notations for hyperbolic functions: 

sh = sinh(sv), ch = cosh(sv), sh = sinh(rv), ch = cosh(rv) . (9) 

The first term in (7) corresponds to a constant A4 component at spatial infinity {A4 w iv^) and gives rise to the 
non-trivial holonomy. One can see that Af^ is periodic in time with period 1 (since we have chosen the temperature 
to be equal to unity). A useful formula for the field strength squared is [4] 

TvF^^F^^ = d^dHogij. (10) 

In the situation when the separation between dyons ri2 is large compared to both their core sizes i (M) and = (L), 
the caloron field can be approximated by the sum of individual BPS dyons, see Figs. 1,3 (left) and Fig. 4. We give 
below the field inside the cores and far away from both cores. 



A. Inside dyon cores 

In the vicinity of the L dyon center zi and far away from the M dyon (sv ^ 1) the field becomes that of the L 
dyon. It is instructive to write it in spherical coordinates centered at Zi. In the 'stringy' gauge [9] in which the A4 
component is constant and diagonal at spatial infinity, the L dyon field is 

= ^ (^^ + 27r - vcoth(vr)^ , A^ = 0, 

._— sin(27ra;4 — </>) n + cos(27ra;4 — (/))t2 

= 2sinh(vr) ' 

, r cos(27ra;4 — 6) tt + sin(27ra;4 — (i)T2 tan(^/2) , , 

Aa = tv — \; . , , — 7^ — ^-^-zrs . (11) 

2sinh(vr) 2r ^ ' 



FIG. 3: The action density of tlie KvBLL caloron as function of z, x at fixed t — y — 0. At large separations ri2 the caloron 
is a superposition of two BPS dyon solutions (left, ri2 = 1.5/r). At small separations they merge (right, ri2 — 0.6/T). The 
caloron parameters are the same as in Fig. 1. 



Here Ag , for example, is the projection of A onto the direction fig — (cos 9 cos cj), cos 6 sin (j), — sin 6) . The 4> component 
has a string singularity along the z axis going in the positive direction. Notice that inside the core region (vr < 1) 
the field is time-dependent, ahhough the action density is static. At large distances from the L dyon center, i.e. far 
outside the core one neglects exponentially small terms O {e~^^) and the surviving components are 



r — *oo 1^^112 t/ 

^ r 2 

corresponding to the radial electric and magnetic field components 

1 ir^ 
~2 



E^^-By-;~y-. (12) 



This Coulomb-type behavior of both the electric and magnetic fields prompts the name 'dyon'. 

Similarly, in the vicinity of the M dyon and far away from the L dyon (rv 3> 1) the field becomes that of the M 
dyon, which we write in spherical coordinates centered at £2: 



a/ _ ^'^3 / ^ 1\ ,M sin0 Ti - COS0 r2 



vcoth(vs)-- , = 



2 \ s / ' 2i sinh(vs) 



iM_n _ cos(/) Ti sin(/) r2 . tan(6'/2) 



whose asymptotics is 



=«' 2.s:nh(v.) (1^^ 



.If (^v-i 

ijf = Bf (14) 

We see that in both cases the L,M fields become Abelian at large distances, corresponding to (electric, magnetic) 
charges (— ,— ) and (+,+), respectively. The corrections to the fields (11,13) are hence of the order of l/ri2 arising 
from the presence of the other dyon. 



B. Far away from dyon cores 

Far away from both dyon cores (rv ^1, sv ^ 1; note that it does not necessarily imply large separations - the 
dyons may even be overlapping) one can neglect both types of exponentially small terms, O {e~"\ and O [eT'^^). 
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With exponential accuracy the function x in eq.(8) is zero, and the KvBLL field (7) becomes Abelian [4]: 

^r = *y (W + C^^'ln^,,), (15) 

where <l?as is the function $ of eq.(8) evaluated with the exponential precision: 

r + s + ri2 S-S3 ^ ^ 

*as = — ; = if ri2 = ri2e3. (16) 

r + s — ri2 r — rz 

It is interesting that, despite being Abelian, the asymptotic field (15) retains its self-duahty. This is because the S'* 
color component of the electric field is 

Ef = diAl = didaln 

while the magnetic field is 

Bf = tijkdjAl = dA In - Sad^ In $as, 

where the last term is zero, except on the line connecting the dyon centers where it is singular; however, this singularity 
is an artifact of the exponential approximation used. Explicit evaluation of eq.(15) gives the following nonzero 
components of the field far away from both dyon centers: 

^r = ^fv+i-iV (17) 



2 V r s 



Z73 / 1 _^ 1\ / (ri2 - r + .s)(ri2 + r - s) 
2 \r s J y {ri2 + r + s){r + s - ri2) 



In particular, far away from both dyons, A4 is the Coulomb field of two opposite charges. 



III. THE SCHEME FOR COMPUTING DET(-D^) 

As explained in section I, to find the quantum weight of the KvBLL caloron, one needs to calculate the small 
oscillation determinant, Det(— _D^), where = 9^ + A^ and is the caloron field (7). Instead of computing the 
determinant directly, we first evaluate its derivative with respect to the holonomy v, and then integrate the derivative 
using the known determinant at v = [2] as a boundary condition. 

If the background field A^ depends on some parameter P, a general formula for the derivative of the determinant 
with respect to such parameter is 



aiogDet(-£»2[A]) 



= -J d^xTr{drA^J^), (19) 



&P 

where is the vacuum current in the external background, determined by the Green function: 

Jf= {6^j'^d,-6^Xdy+A^%+A'''S^,)g'^''{x, y)l^^ or simply = D^G + QD^. (20) 
Here Q is the Green function or the propagator of spin-0, isospin-1 particle in the given background A^ defined by 

-DlG{x,y)=S^'\x-y) (21) 

and, in the case of nonzero temperatures, being periodic in time, meaning that 

00 

G{x,y)= ^ G{x4,x;y4 + n,y}. (22) 

n— — 00 

Eq.(19) can be easily verified by differentiating the identity logDet(— £)^) = Trlog(— £)^) [35]. The background field 
Afj^ in eq.(19) is taken in the adjoint representation, as is the trace. Hence, if the periodic propagator Q is known, 

eq.(19) becomes a powerful tool for computing quantum determinants. Specifically, we take P = v as the parameter 
for differentiating the determinant, and there is no problem in finding d^Ai^ for the caloron field (7). 
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The Green functions in self-dual backgrounds are generally known [25, 26] and are built in terms of the Atiah- 
Drinfeld-Hitchin-Manin (ADHM) construction [28] for the given self-dual field. A subtlety appearing at nonzero 
temperatures is that the Green function is defined by eq.(21) in the Euclidean IR^ space where the topological charge 
is infinite because of the infinite number of repeated strips in the compactified time direction, whereas one actually 
needs an explicitly periodic propagator (22). To overcome this nuisance, Nahm [25] suggested to pass on to the 
Fourier transforms of the infinite-range subscripts in the ADHM construction. We perform this program explicitly in 
Appendix A, first for the single dyon field and then for the KvBLL caloron. In this way, we get the finite-dimensional 
ADHM construction both for the dyon and the caloron, with very simple periodicity properties. Using it, we construct 
explicitly periodic propagators in Appendix B, also first for the dyon and then for the caloron case. For the KvBLL 
caloron it was not known previously. Using the obtained periodic propagators, in Appendix C we calculate the exact 
vacuum current (20) for the dyon, and in Appendix D we evaluate the vacuum current in the caloron background, 
with the help of the rcgularization carried out in Appendix E. 

Although there is no principle difficulty in doing all calculations exactly for the whole caloron moduli space, at some 
point we loose the capacity of performing analytical calculations for the simple reason that expressions become too 
long, and so far we have not been able to put them into a manageable form in a general case. Therefore, we have to 
adopt a more subtle attitude. First of all we restrict ourselves to the part of the moduli space corresponding to large 
separations between dyons (ri2 ^ 1). Physically, it seems to be the most interesting case, see section L Furthermore, 
at the first stage we take ri2V, ri2V ^ 1, meaning that the dyons are well separated and do not overlap since the 
separation is then much bigger than the core sizes, see Figs. 1,3 (left). In this case, the vacuum current (20) becomes 
that of single dyons inside the spheres of some radius R surrounding the dyon centers, such that k ^ i? <c ri2, and 
outside these spheres it can be computed analytically with exponential precision, in correspondence with subsection 
II. B, see Fig. 4. Adding up the contributions of the regions near two dyons and of the far-away region, we get 
dDet{—D^)/dY for well-separated dyons. Integrating it over v we obtain the determinant itself up to a constant and 
possible l/ri2 terms. 

This is already an interesting result by itself, however, we would like to compute the constant, which can be 
done by matching our calculation with that for the trivial caloron at v = 0. It means that we have to extend the 
domain of applicability to ri2V = C(l) (or ri2V = C(l)) implying overlapping dyons, presented in Figs. 1,3 (right). 
To make this extension, we 'guess' the analytical expression which would interpolate between ri2V ^ 1 where the 
determinant is already computed and ri2V <C 1 where matching with the Gross-Pisarski-Yaffe (GPY) calculation [2] 
can be performed. At this point it becomes very helpful that we possess the exact vacuum current for the caloron, 
which, although too long to be put on paper, is nevertheless affordable for numerical evaluation (and can be provided 
on request). We check our analytical 'guess' to the accuracy better than one millionth. In this way we obtain the 
determinant up to an overall constant factor for any v, v with the only restriction that ri2 ^ 1 . This constant factor 
is then read off from the GPY calculation [2]. 

Finally, we compute the l/ri2 and logri2/r'i2 corrections in the Det(— I?^), which turn out to be quite non-trivial. 
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IV. DET(-D^) FOR WELL SEPARATED DYONS 

The L,M dyon cores have the sizes = and i, respectively, and in this section we consider the case of well-separated 
dyons, meaning that the distance between the two centers is much greater than the core sizes, ri2 ^ ^, This 
situation is depicted in Figs. 1,3 (left). The two dyons are static in time, so that d logDet(— _D^)/9v (19) becomes 
an integral over 3d space, times 1/T set to unity. We divide the 3d volume into to three regions (Fig. 4): i) a ball of 
radius R centered at the center of the M dyon. ii) a ball of radius R centered at the L dyon, iii) the rest of the space, 
with two balls removed. The radius R is chosen such that it is much larger than the dyon cores but much less than 
the separation: ri2 ^ ^, =■ Summing up the contributions from the three regions of space, we are satisfied to 
observe that the result does not depend on the intermediate radius R. 



A. Det(— D"^) for a single dyon 



In region i) the KvBLL caloron field can be approximated by the M dyon field (13), and the vacuum current by that 
inside a single dyon, both with the C(l/ri2) accuracy. We make a more precise calculation, including the 0{l/ri2) 
terms, in section V. The single-dyon vacuum current is calculated in Appendix C. Adding up the three parts of the 
vacuum current denoted there as J^'"''™ we obtain the full isospin-1 vacuum current (in the 'stringy' gauge) 



Jr = 0, 



Jtj) = — 

Je = - 

Ji = -in 

~acb 



iv(sin(^)r2 + cos(0)Ti)(l - svcoth(sv))2 



247^2^2 ginh(sv) ' 
zv(sin(0)ri - cos(0)r2)(l - svcoth(sv))2 
247r^s^ sinh(sv) ' 
(1 — svcoth(sv))^ 1 — svcoth(sv) coth(sv)(l — svcoth(sv)) 



+ 



3s 



+ ■ 



27rs2 



where (Tc)° = ie"''* are the isospin-1 generators. We contract (23) with dA^/dY from eq.(13) according to eq.(19). 
After taking the matrix trace, the integrand in eq.(19) becomes spherically symmetric: 



Tr[avA^J^] 



2 

3s 



coth(sv) — sv + 



sv (svcoth(sv) — 2) 



sinh (sv) 

(svc:oth(.s-v) - I)'"' (sinh(2.sv) - 3.w) (sv cot,h(.sv) - 1)" (siuh(2.svj - 2»\') 



67r^ sinh^ (sv)s^ 



27r sinh^(sv) tanh(sv)s^ 



(23) 



It has to be integrated over the spherical box of radius R. Fortunately, we are able to perform the integration 
analytically. The result for the M dyon is 



9 log Pet (-D^[M dyon]) 



^ r , o 24(7B-log7r 



53 + Itt^ 



IStt 



+ 



1 



+ 



^ P' (v) - 27ri?2 P" (v) + 2TrR P'" (v) - — log(i?v) . 
3 37r 



(24) 



As we see, it is badly infrared divergent, as it depends on the box radius R. Here P{q) is the potential energy (see 
eq.(2)) 



P{q) = 



12 Vtt / Vtt 



P'{q) = ^q{7r-q){27r-q), 



P"{q) = :^{3q^-67Tq + 2n^), p'" (q) = (q - , P^^ {q) = 

The IR-divergent terms arise from the asymptotics of the integrand. Neglecting exponentially small terms e~*'^ in 
eq.(23) we have 



-Tvid^A^J^] 



(1 - SV)3 (1 - SV)2 1 - sv 



127r2s3 



+ 



+ ■ 



6s 



P' (v - 1) = P'(v) - P"(v)i + Ip-'M^ - IP",,)^. 
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Integrating it over the sphere of radius R one gets the IR-divergent terms (the second Hnc in cq.(24)). 

The fact that the IR-divergent part of dDet{—D^)/dv is directly related to the potential energy P(^4) is not 
accidental. At large distances the field of the dyon becomes a slowly varying Coulomb field, sec eq.(14). Therefore, 
the determinant can be generically expanded in the covariant derivatives of the background field [12, 13] with the 
potential energy P(^4) being its leading zero-derivative term. The nontrivial fact, however, is that with exponential 
precision the vacuum current is related to the variation of solely the leading term in the covariant derivative expansion 
of the eflFective action with no contribution from any of the subleading terms. This is a specific property of self-dual 
fields, and we observe it also in the following subsection. 



B. Contribution from the far-away region 



We now compute the contribution to dDet{— D"^) / dv from the region of space far away from both dyon centers. 
With exponential accuracy (meaning neglecting terms of the order of e~'"^ and e"*^) the KvBLL caloron field is given 
by eqs.(17,18), and only the A4 component depends (trivially) on v. The caloron vacuum current with the same 
exponential accuracy is calculated in Appendix D. Combining the results given by eqs.(D8,D9) and eqs.(Dll,D12) we 
see that J„ = and for J4 we have 



J4 



El 
2 



-I- 



( 








37r2 ' 


,r s) 


/I 1^ 




- 




1 - 1677 


\r s y 





i i_i 

n \r s 



, 647r o 



2 /I 1 



3 \r 



4 /I 1 
TT \ r s 



1 1 

r s 



»7r 

T 



(25) 



We remind the reader that r, s are distances from M,L dyon centers 2*1,2 and that u) = v/(47r). It is interesting that 
the separation ri2 = \zi — Z2\ does not appear explicitly in the current. Moreover, it can be again written through 
the potential energy ^(^4): 



J4 = ^iTs P'{q)\g^y+i/r-l/s ■ 



Therefore, in the far-away region one obtains 



-Tr[d^A^J^,]=P' (v 



1 1 
r s 



(26) 



(27) 



We have now to integrate eq.(27) over the whole 3d space with two spheres of radius R surrounding the dyon centers 

removed: 



(28) 



- j d^xTrld^A^J^,] = Jd^xP' (y+^-^^ 

= P'(v) / d^x + P"(v) jd^x (i - i) + iF'"(v) jd^x + ^^^^(v) jd^x [I - I 

The first integral in eq.(28) is the 3d volume V, minus the volume of two spheres, V — 2^R^. The second integral is 
zero by symmetry between the two centers, and so is the last one. The only non-trivial integral is 



47rri2 - 87ri? + 0(^1. 



Therefore, the contribution from the region far from both dyon centers is 



d log Betj-D'^ 



far 



An _ 



1 



P'(v) V - + -P"'(v) (4^ri2 - 87ri?) . 



(29) 



(30) 
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C. Combining all three regions 

We now add up the contributions to 9 log Det(— -D^)/9v from the regions surrounding the two dyons and from 
the far-away region. Since the contribution of the L dyon is the same as that of the M dyon with the replacement 

V ^ V and since d/dv = —d/dv, when adding up contributions of L,M core regions we have to antisymmetrize in 

V -(-^ V. It should be noted that P(v) and P"{v) are symmetric under this interchange, while P'(v) and P"'{v) are 
antisymmetric. Therefore, the combined contribution of both cores is, from eq.(24), 



d log Det(-£'2) 



dv 



2P'(v)^i?3 + 1 STTi? + 1 - 1 - f In (^) . (31) 

3 2 V V 37r Vv/ 



Adding it up with the contribution from the far-away region, eq.(30), we obtain the final result which is independent 
of the intermediate radius R used to separate the regions: 

d log Det(-D^) ^ ^ ^ ^,„^^^ ^^^^^ ^ 1 _ 1 4 ^ , ^3^^ 



9v V V 37r Vv 

This equation can be easily integrated over v up to a constant which in fact can be a function of the separation ri2: 

log Det(-£>2) = p(v) V + P" (v) 27rri2 + (l - log(v) + " £) ^°s(v) + /(n2) • (33) 

Since in the above calculation of the determinant for well-separated dyons we have neglected the Coulomb field of one 
dyon inside the core region of the other, we expect that the unknown fmiction f{r\2) = 0{l/r\2) + c, where c is the 
true integration constant. Our next aim will be to find it. The 0{l/ri2) corrections will be found later. 



V. MATCHING WITH THE DETERMINANT WITH TRIVIAL HOLONOMY 

To find the integration constant, one needs to know the value of the determinant at v = (or v = 0) where the 
KvBLL caloron with non-trivial holonomy reduces to the Harrington-Shepard caloron with a trivial one and for which 
the determinant has been computed by GPY [2]. Before we match our determinant at v^^O with that at v = let us 
recall the GPY result. 



A. Det(-D^) at V = 

The V = periodic instanton is traditionally parameterized by the instanton size p. It is known [2, 24] that the 

periodic instanton can be viewed as a mix of two BPS monopoles one of which has an infinite size. It becomes 
especially clear in the KvBLL construction [4, 5] where the size of one of the dyons becomes infinite as v ^ 0, see 
section II. Despite one dyon being infinit ely large, one can still continue to parametrize a caloron by the distance ri2 
between dyon centers, with p = ■y/ri2/7r. Since our determinant (33) is given in terms of ri2 we have first of all to 
rewrite the GPY determinant in terms of ri2, too. Actually, GPY have interpolated the determinant in the whole 
range of p (hence ri2) but we shall be interested only in the limit ri2 » 1. In this range the GPY result reads: 

logDet(-l)2)|^^^^^^^ = logDet(-D2)|^^^_^^^ + l7rri2-^logri2 + Co + o(^), (34) 



Co = x 



8 87ij 2 7r^ 4 1og7r 

9 3 Tf^ 3 



We have made here a small improvement as compared to ref. [2], namely i) we have checked that the correction is 
of the order of 1/ r\2 basing on an intermediate exact formula, ii) we have also managed to get an exact analytical 
expression for the constant. 

The zero-temperature determinant is that for the standard BPST instanton [19, 27]: 

logDet(-D2)|^^^_^^^ = ^log^+ilog(^)+a(l), (35) 

2^E 16 , log 2 , 21og(27r) 4C'(2) 
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where it is implied that the determinant is regularized by the Pauli-Villars method and ji is the Pauli-Villars mass, 
see section VI. A. Combining eqs. (34,35) one obtains 

logDet[-£»2]|^^o^ = ^7rri2 - logns + | log/. + c^+O (^-^^ , (37) 

where 

ci = log 2 + ^ log TT - 3 - 2 7B - ^ - ^^4^ = 0.206602292859 . (38) 

We notice that P{v) and P"{v) ^ | at v — !■ 0, therefore the first two terms in cq.(33) become exactly equal to 
the first term in eq.(37). At the same time, the last two terms in eq.(33) become log v— | log(27r) +c which is formally 
singular at v ^ and does not match the — logri2 in eq.(37). The reason is, eq.(33) has been derived assuming 
^12 ^ 7; 1^ ™d oiic cannot take the limit v ^ in that expression without taking simultaneously ri2 ^ oo. In order 
to match the determinant at v = one needs to extend eq.(33) to arbitrary values of vri2. As we shall see, it will be 
important for the matching that logri2 has the coefficient —1. 



B. Extending the result to arbitrciry values of vri2 



Let us take a fixed but large value of the dyon separation ri2 » 1 such that both eq.(33) and eq.(37) are valid. 
Actually, our aim is to integrate the exact expression for the derivative of the determinant 



avlogDet(-L>2) = j p{x)d*x, p{x) = -Tr[d^Ai,J''{x,x\A)], 



(39) 



from v = where the determinant is given by eq.(37), to some small value of v ^ 1 (but such that vri2 ^ 1) where 
eq.(33) becomes valid. We shall parametrize this v as v = fc/ri2 "C 1 with k ^ 1. The result of the integration over 
V must be equal to the difference between the right hand sides of eqs. (33,37). We write it as 



d'^x p{x) = V 



P(A)_p(o) 



+27rri2 



P"(^)-P"(0) 



+log(A;)+c--log(27r)-ci--log/x+e» — 



Notice that logri2 has cancelled in the difference in the r.h.s. We denote 

2 5 
C2 = c - ci - - log/i - - log(27r) . 



r-i2y 
■ (40) 



(41) 



We know that the first two terms in eq.(40) come from far asymptotics. Denoting by p our p with subtracted 
asymptotic terms we have 



d'^x p{x) = logfc + C2 + C ( — 

ri2 



(42) 



In this integration we are in the domain 1/v ^ 1 and ri2 ^ 1 and we can simplify the integrand dropping terms 
which are small in this domain. At this point it will be convenient to restore temporarily the temperature dependence. 
With (3 = 1/T our domain of interest is 1/v ^ /? and ri2 ^ /3. Therefore we are in the small-/? domain and can 
expand p in series with respect to /3: 



P=^Po + ^Pi + 0(/?°). 



(43) 



As we shall see in a moment, only the first two terms are not small in this domain and we need to know only them to 
compute C2. It is a great simplification because po,i do not contain terms proportional to e^^'' since v = 2itT — v oo 
at /3 ^ 0, and what is left is time independent. Moreover, what is left after we neglect exponentially small terms are 
homogeneous functions of r, s, ri2, v and we can turn to the dimensionless variables: 



Po{r,s,ri2,Y) 
Pi(r-,s,ri2,v) 



1 _ / r s V 

— Po — , — , — 
r-12 \ri2 r-12 ri2 

1 ~ / r s V 

— Pi — , — , — 
fi2 V^i2 ri2 ri2 
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We rewrite the l.h.s. of eq.(42) in terms of the new quantities: 



J d^x p{x) = ^ Jd^xpo+ ri2 J d^x 



+ ri2 J d^x pi + 0{(3) , (44) 

where x = x/ru is dimensionless. We see that it is indeed sufficient to take just the first two terms in the expansion 
(43) at /3 — > 0. The integration measure can be written in terms of the dimensionless variables r = r/ri2, s = sjrxi 
as 

(^x = 2-K fdf s ds , (45) 

where f and s arc constrained by the triangle inequalities f+s<l, f+l<s and s+l<f, and we have integrated over 
the azimuth angle. 

We have now to use the exact vacuum current to compute / po,i- First, it turns out that the first integral in eq.(44) 
is zero. This is good news because had it been nonzero, eq.(42) could not be right as its r.h.s. has no dependence on 
ri2 other than possible l/ri2 terms. Second, we have noticed that the second integral in eq.(44) is in fact 



/ 



d^x pi = . (46) 

vri2 + 1 



Unfortunately, we were not able to verify it analytically but we checked numerically that it holds with the precision 
of a few units of 10"'' in the range of vri2 between and 15. Combining eqs. (44,46) we obtain for the l.h.s. of eq.(42) 

^ dv J d^x p{x) = ri2 ^ = \og{k + 1) = log fc + O . 

Therefore, we reproduce the r.h.s. of eq.(42) and in addition find that C2 = 0. 

Eq.(46) is sufficient to extend the result for the determinant (33) valid at vri2 1 to arbitrary values of vri2, 
provided ri2 '> 1 (the extension to arbitrary values of vri2 is obtained by symmetry v ^ v). The final result for the 
determinant to the l/ri2 accuracy is 

logDet[-D2] = VP{v) + 27rF"(v) ri2 + (l " ^) log (^^12 + 1) + " (^^12 + 1) 

+ ^log(/iri2)+ci + |log(27r) + 0( — ) (47) 



3 3 \ri2 

where /z is the UV cutoff and the numerical constant ci is given by eq.(38). This expression is finite at v ^ 0, v ^ 
and coincides with the GPY result (37) in these limits. At vri2 ^ 1 we restore the previous result, eq.(33), but 
now with the integration constant fixed: c = |log/i + | log(27r) + Ci. Eq.(47) is valid for any holonomy, i.e. for 
V, V e [0, 2tt], and the only restriction on its applicability is the condition that the dyon separation is large, ri2 1- 



C. l/ri2 corrections 

Eq.(47) can be expanded in inverse powers of vri2, vri2, which gives l/(vri2), l/(vri2) (and higher) corrections; 
however, there are other l/ri2 corrections which are not accompanied by the 1/v, 1/v factors: the aim of this 
subsection is to find them using the exact vacuum current. 

To this end, we again consider the case ri2 3> -^j = such that one can split the integration over 3d space into three 
regions shown in Fig. 4. In the far-away region one can use the same vacuum current (25) as it has an exponential 
precision with respect to the distances to both dyons. In the core regions, however, it is now insufficient to neglect 
completely the field of the other dyon, as we did in section IV looking for the leading order. Since we are now after 
the l/ri2 corrections, we have to use the exact field and the exact vacuum current of the caloron but we can neglect 
the exponentially small terms in their separation. 

Another modification with respect to section IV is that we find it more useful this time to choose ri2 as the 
parameter V in eq.(19). We shall compute the l/rj2 terms in dDct{—D'^)/dri2 and then restore the determinant itself 
since the limit of ri2 ^ 00 is already known. Let us define how the KvBLL field depends on ri2. As seen from eq.(7) 
the KvBLL field is a function of r, s, v, ri2 only. We define the change in the separation ri2 — » ri2 + dri2 as the 
symmetric displacement of each monopole center by ±c?ri2/2. It corresponds to 

dr rl2+r'^ — s'^ ds + s"^ — r'^ 

dri2 4ri2r ' dri2 4ri2S 
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We shall use the definition (48) to compute the derivative of the caloron field (7) with respect to ri2. 

Let us start from the M-monopolc core region. To get the l/ri2 correction to the determinant we need to compute 
its derivative in the 1 /r^2 order and expand correspondingly the caloron field and the vacuum current to this order. 
Wherever the distance r from the far-away L dyon appears in the equations, we replace it by r = (rf2 + 2sri2 cos 9 + 
s^)!/^ where s is the distance from the M-dyon and 9 is the polar angle seen from the M-dyon center. Expanding in 
inverse powers of ri2 we get coefficients that are functions of s, cos 9. One can easily integrate over 9 as the integration 
measure in spherical coordinates is 27rs^ds dcoii9. We leave out the intermediate equations and give only the end 
result for the integrand in eq.(19). After integrating over cos^ we obtain the following contribution from the core 
region of the M monopole: 



SDetC-f^) 



dr 



12 



M dyon core 



= \- [ /w^2 l&TT S^ds + O f 



(49) 



where Ii/rf^ reads 



coth(sv) coth(sv) sv coth(sv)csch(sv) s v^ (2 + cosh(2sv)) csch(sv) 
" 127r2s3 9s 36 72 



v (-61 + 3 cosh(2sv) + 4 cosh(4sv)) coth(sv)csch(sv)'^ v^ (37 + 23 cosh(2sv) + 4 cosh(4sv)) coth(sv)csch(sv)^ 



967r.s 1927r2s 
s^v^ (4 + cosh(2sv)) coth(sv)^csch(sv)'* sv^ (54cosh(sv) + 17cosh(3sv) + cosh(5sv)) cscli(,sv)^ 

487r 38477 
sv'* (-406 cosh(sv) - 81 cosh(3sv) + 7 cosh(5sv)) csch(sv)'^ sV^ (-24 - 33 cosh(2sv) + cosh(6sv)) csch(sv)* 



23047r2 



^ 6 + 13csch(sv)^ ^ —10 — 5csch(sv)^ + 9csch(sv)'* ^ 40 + 65csch(sv)^ + 19csch(sv) 



15367r2 

4 



+v 



72 ■ 487rs2 ' ' 1927r2s2 

2 -6 - 20csch(sv)^ + 9csch(sv)'* + 27csch(sv)^ 3 24 + 98csch(sv)^ + 285csch(sv)'* + 234csch(sv)^ 



487r 



+ V'' 



5767r2 



(50) 



Fortunately we are able to integrate this function analytically: 



L 



, r 1. 2, 1 7r2 + 367i; + 69 , 2v(v2 - 37rv + 27r2) ^3 , 4(6v7r - 27r2 - 3v2) , 10(v - tt) „ 41og(i?v/7r) 
'0 A/.?,167rs ds = — + R + R + 3^ ^ • 

(51) 

For the L monopole core contribution one has to replace v by v. Adding together contributions from L,M monopole 
cores we have 



aDet(-£>2) 



dr 



12 



1 



'12 



1 1 _ ,^- + 367. + 69 _ ^ _ ^^2)^2 _ 41og(^!v^ 

V V 277r 97r^ ^ 3w 



(52) 



Now let us turn to the far-away region. Recalling eq.(26) we realize that the contribution of this region is determined 
by the potential energy: 



aDet(-£)2) 



dr 



12 



far 



d^xdr.^P v + 



1 



-P"(v) d^xd, 



1 1 

r s 



1 



1 1 

r s 



The integration region is the 3d volume with two balls of radius R removed. We use 



1 1 

r s 



2n 
Wo. 



4SI„.(!|) 



97r^ + 8 



(53) 



Adding up all three contributions we see that the region separation radius R gets cancelled (as it should), and we get 



dr,, logDet(-D2) = 27rP"(v) + — 



12 



4 , 



' 12 



_ 1 _ 1 _^ 50 _^ 87_E 237r 
V V 97r 377 54 



(54) 



which can be easily integrated, with the result 



logDet(-£»2) = 27rP"(v) na + — 

ri2 



1 _^ 1 ^ 237r 87E 74 4 
V V 54 37r 97r 37r 



' 12 



+ -c + o[-^J (55) 
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where c is the integration constant that does not depend on ri2. Comparing eq.(55) with eq.(32) at ri2 — > oo we 
conclude that 

2 Stt — 4v Stt — 4v 5 

c = VP{v) + - log/i + logv + logv + - log(27r) + ci (56) 

and ci is given in eq.(38). We note that the leading correction, logri2/ri2, arises from the far-away region and is related 

to the potential energy, similar to the loading ri2 term. The terms proportional to i and = can be extracted from 
expanding eq.(47) (which is an additional independent check of eq.(46)). In fact, eqs.(47) and (55) are complementary: 
eq.(47) sums up all powers of but misses ^°f^^^ and ^ terms, whereas eq.(55) collects all terms of that 

order but misses higher powers of ^— , 

o t- ri2V' ri2V 



VI. QUANTUM WEIGHT OF THE KVBLL CALORON 

A. Quantum weight of a Euclidean pseudoparticle: generalities 

If a field configuration is a solution of the Yang-Mills Euclidean equation of motion, D^Fn^ = 0, its quantum 
weight is the contribution of the saddle point to the partition function 

2 = I DA, eM-S[A]), S[A] = ^ j d'x F^^F^,. (57) 

The general field over which one integrates in eq.(57) can be written as 

Af, = Af, + (58) 

where A, is the classical solution corresponding to the local minimum of the action and a, is the presumably small 
quantum oscillation about the solution. One expands the action around the minimum, 



S[A]=S[A]-^ jd'xalDf{A)Fl^{A) + ^ ld^xa;w;t{A)ai + 0{a' 



), (59) 



where the linear term is in fact absent since A satisfies the equation of motion, and the quadratic form is 

W^tiA) = -D\Ar\. + {DM''''{A)-2r>>F^^,{A), (60) 
Df{A) = df.S"^ + f'^^AI. (61) 

We have written the covariant derivative in the adjoint representation; the relation with the fundamental representation 
is given by a, = —ia'^t"', Tv{t°'t^) = ii5"^ and similarly for F,^, etc. The 1-loop approximation to the quantum weight 
corresponds to evaluating eq.(57) in the Gaussian approximation in a^, hence C(a^) terms in eq.(59) have been 
neglected. 

The quadratic form (60) is highly degenerate since any fluctuation of the type aj^ = D1^[A)A!'{x) corresponding to 
an infinitesimal gauge transformation of the saddle-point field A, nullifies it. Therefore, one has to impose a gauge- 
fixing condition on a,. The conventional choice is the background Lorenz gauge D1^{A)a^^ = 0: with this condition 
imposed the operator W simplifies as the second term in eq.(60) can be dropped. Fixing this gauge, however, brings 
in the Faddeev-Popov ghost determinant Det(— £'^(A)). 

To define the path integral, one decomposes the fluctuation field in the complete set of the eigenfunctions of the 
quadratic form, 

<(^) = EC"<"(^)' Kl'i'ln = ^ni'ln. « „ = 0, (62) 

n 

and implies that the path integral is understood as the integral over Fourier coefficients in the decomposition: 

DA,{x)=Xl^, (63) 



The quadratic form (60) has a flnite number of zero modes related to the moduli space of the solution. Let the number 
of zero modes be p (for a self-dual solution with topological charge one p = AN for the SU{N) gauge group [26]). Let 



16 



^i, i = be the set of collective coordinates characterizing the classical solution, of which the action S[A] is in 

fact independent. The zero modes are 

r,i{x) = - Df{A)k\{x) (64) 

where the second term is subtracted in order for the zero modes to satisfy the background Lorenz condition, D^^^'tp'^ ^ = 
0. The px p metric tensor 



gij= Jd^xr^iij^j (65) 



defines the metric of the moduli space. Its determinant is actually the Jacobian for passing from integration over 
zero-mode Fourier coefiicients Cj, i = l...p, in eq.(63) to the integration over the collective coordinates ^i, i = l...p: 

n7r=-^n^^^(4=)'' j=v^- (66) 

i=l * 1=1 ^ * ^ 

Finally, one has to normalize and regularize the ghost determinant Det(— -D^) and the Gaussian integral of the 
quadratic form. One usually normalizes the contribution of a pscudoparticle to the partition function by dividing 
it by the free (i.e. zero background field) determinants, and regularizes it by dividing by the determinants of the 
—D"^ and operators shifted by the Pauli-Villars mass fi [19, 27]. It means that Det{—D^) is replaced by the 
'quadrupole' combination 

Bet(-D^) Det(-ZP^) Det(-9^ + M^) 
Uet( V Det(-a2) Det(-£>2 + ^2) 

and similarly for the determinant of the quadratic form, 

n..'m/ ^ Det^(;y^.) Det{-d%^ + 1,^) 

Bet (W^,.)n,r = T. \ r. 4-' fur . — 2\~' (68) 

where the prime indicates that only the product of nonzero eigenvalues is taken. In the integration over Pauli Villars 
fields, the zero eigenvalues are shifted by /U^. Hence the integration over the zero-mode Fourier coefficients in the 
Pauli-Villars fields produces the factor 



/i 



(69) 



which has to be taken in the minus first power. Finally, one obtains the following normalized and regularized expression 
for the 1-loop quantum weight of a Euclidean pscudoparticle: 



e-^W i^^X J (Det'(W^^,)„,,) ^ Det(-£)2)n,r. (70) 



.gV27r 



If the saddle-point field is (anti)self-dual there is a remarkable relation between the two determinants [20]: 
Det'(W^^)n,r = Det^(— D^)n, r which is satisfied if the background field is decaying fast enough at infinity and the 
Hilbert space of the eigenfunctions of the two operators is well defined. This is the case of the KvBLL caloron but not 
the case of a single BPS dyon having a Coulomb asymptotics. To define the dyon weight properly, one would need to 
consider it in a spherical box, which would violate most of the statements in this subsection. For this reason we prefer 
to consider the well-defined quantum weight of the KvBLL caloron in which case the product of two determinants in 
eq.(70) becomes just T>et~^{-D'^). 



B. KvBLL caloron moduli space 



The KvBLL moduli space has been studied in the original papers [4, 5]; in particular in ref. [4] the metric tensor 
Qij (65) has been explicitly computed. We briefiy review these results and adjust them to our needs. 
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The KvBLL classical solution has 8 parameters for the SU{2) gauge group. These are the four centcr-of mass 
positions and the four quaternionic variables ( = pU corresponding to the constituent monopoles relative position 
in space and one global gauge transformation, see Appendix A. 3. The moduli space of the KvBLL caloron is a product 
of the base manifold IR"^ x parameterized by the z G IR"^ and 24 e [0, 1], and the non-trivial part of the moduli 
space parameterized by the quaternion Q. It should be noted that the change C, — corresponding to the center of 
the SU{2), leaves Afj_{x) invariant, such that one has to mod-out this symmetry. 

The 8 zero modes ^ (64) satisfying the background Lorenz condition have been explicitly found in ref. [4] . If one 
parametrizes the unitary matrix through Euler angles, 



the metric is [4] 



where 



{7 = g-iT^g^(f-^e)^g-^v:l^ 0<T<47r, < < 27r, O<0<n, (71) 
rfs^ = {2'Kf[2dz^dz^ + (1 + Stt^wwp^) {Mp^ + p^cfn) + p'^{l + ^n'^uiujp'^y^dYll] (72) 
d'^n = sm^e d>f^ + rfSs = dT + cos 6 dip. (73) 



The first part describes the flat metric of the base manifold x S^, the remainder forms the non-trivial part of 
the metric. The variables are inside the ranges p € [0, 00), 9 € [0,7r), </> € [0, 27r), T S [0, 47r)/.^2 = [0, 27r) for the 
non-trivial part, and Z4 G [0, 1], Zi £JR for translational modes. 

The collective coordinate Jacobian is immediately found from eq.(72): 



J = ^Jdet{gij) = 8 (27r)^ p^ (1 + Sn^p^uu)) sin . (74) 



The factor sin^ is needed to organize the orientation SO{3) Haar measure normalized to unity, 



d^O=— J dT j dipj d0 sm0 = 1, (75) 
and the KvBLL measure written in terms of the caloron center, size and orientation becomes 

j d^z j dzij d^O j dpp^ (1 + Sn^Loujp^) 16 (27r)^°. (76) 

This must be multiplied by the factors Q^/^ ) and exp(— SfA]) = exp ^— according to eq.(70). As the result, 
the KvBLL measure is 

J d'z J dz,j d'0 1^ (1 + Sn'ivQp') (^)' (77) 

When the holonomy is trivial (w = or a) = 0) it becomes the well-known measure of the BPST instanton [27] or that 
of the Harrington-Shepard caloron [2] . The difference between the two is that in the first case one integrates over any 
Z4 whereas in the second case the Z4 integration is restricted to Z4 G [0, Eq.(77) would have been the full result in 
the M=l supersymmetric theory where the determinant over nonzero modes is cancelled by the gluino determinant. 
In that case one would need to add the integral over Grassmann variables corresponding to the gluino zero modes. 



C. Combining the Jacobian and the determinant over nonzero modes 

According to the general eq.(70), we have now to multiply eq.(77) by the (regularized and normalized) determinant 
over nonzero modes, which has been calculated in eq.(47). First of all, we notice that Det~^(— D^) brings in an 
additional UV divergent factor /U~3. In combination with the classical action and the factor p,^ coming from zero 
modes, it produces 

22 22 

p~ e = A~ (78) 
where A is the scale parameter obtained here through the 'transmutation of dimensions'. 
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We notice further that Det~"'^(— D^) is independent of the SU{2) orientation O and of Z4. Therefore, we integrate 
over these variables, which gives unity. Next, we introduce the centers of the constituent BPS dyons £1^2 such that 
\zi — Z2\ = ri2 = Trp^ and write 



Zl + Z2 



d-^{zi - Z2) = 47r / (rzr{2drY2 = Stt^ / d^zdpr 



(79) 



Therefore, integration over d'^zdp in eq.(77) can be traded for integrating over the dyon positions in space, Zi^2- Lastly, 
we restore the temperature from dimensional considerations and obtain our final result for the 1-loop quantum weight 
of the KvBLL caloron, written in terms of the coordinates of the dyon centers: 



■'KvBLL 



-I 



d^zi d^Z2 C 



Stt^ 



AttT 



Tri2 



2tt- 



vv 



ri2 ) (vri2 + 1) ^^^■^ ^ (vri2 + 1) ^""^ ^ 



X exp [-VP{y) - 277 ri2 P"(v) 



where 



C = exp 



8 16 7e 



2^ 
"27 



4C'(2) 



= 1.031419972084 



and P(v) is the potential energy 
1 



P(v) 



127r2r 



2—2 
-V V , 



P"(v) 



7r2r 



ttT 1 



1 

71 



2-kT- 



(80) 



(81) 



(82) 



We have collected the factor Aire ^'^T/A because it is the natural argument of the running coupling constant at 
nonzero temperatures [12, 30]. Here A is the scale parameter in the Pauli-Villars regularization scheme that we have 

used. It is related to scale parameters in other schemes: Apy = e^A,^ = 40.66 • exp ^— ■'^lat [31]- The factor 

is not renormalized at the one-loop level: it starts to 'run' at the 2-loop level, see below. 

The KvBLL caloron weight (80) has been derived assuming the separation between constituent dyons is large in 

temperature units (ri2 S> j^) but the holonomy is arbitrary: |TrL <E [—1,1] corresponding to v,v € [0, 27rT]. It 
means that eq.(80) is valid not only for well-separated but also for overlapping dyons. 



D. The limit of large dyon separation 
In the limit when the separation of dyons is larger than their core sizes, ri2 '> ^ , = , the caloron weight simplifies to 



iKvBLL = exp 

X exp 



3 ^ ' 



d^zi d^Z2 (27r)3 C 



8^2 



V 47rr ) 



j,f-(l_j.)f(i-) 



-27rri2T [--Av{l-v) 



(83) 



where we have introduced the dimensionlcss quantity y = € [0, Ij. In subsection V.C we have calculated the 

correction to the determinant, see eq.(55). Another correction arises from the Jacobian (74) which cancels the ^, | 
terms in eq.(55). As a result, we get the following correction factor to eq.(83) 



exp 



1 



ri2r V37r 



■log [Kl-i^)(2ri2T)2] +ci/,,, 



O 



{n2Ty 



Cl/ri2 



74 SjE 237r 



97r 37r 



54 



= 1.946 . (84) 



One can define the interaction potential between L, M dyons as 

VLM(ri2) = ri2 277 Q - 4 y{l - i')) - ^ (J^ log [^^(1 - t^){2n2Tf] + Ci/,,,) + ^ 



(85) 



This interaction is a purely quantum effect: classically L, M dyons do not interact at all as the KvBLL caloron of 
which they are constituents has the same classical action for all X, M separations. Curiously, the interaction potential 
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has the famihar "hnear + Coulomb" form. Both terms depend seriously on the holonomy: the Polyakov line at 
spatial infinity is ^TrL = cos(7ri'). In the range 0.787597 < ^iTrLj < 1 dyons experience asymptotically a constant 
attraction force; in the complementary range 5|Tri| < 0.787597 it is repulsive. It should be noted that in its domain 
of applicability ri2 > second term in eq.(85) is a small correction as compared to the linear rising 

(or linear falling) interaction. 



E. 2-loop improvement of the result 

The factor in eq.(80) is the bare coupling which is rcnormalizcd only at the 2-loop level. In the case of 

the zero-temperature instanton one can unambiguously determine the 2-loop instanton weight without explicit 2-loop 
calculations from the requirement that it should be invariant under the simultaneous change of the UV cutoff and of 
the bare coupling given at that cutoff, such that the scale parameter 

remains invariant. The result [32] is that one has to replace the combination of the bare coupling constants 
where 

/3(r)=6ilnJ, = /3(r) + ^ In (88) 

and r is the scale of the pseudoparticle, which is l/p in the instanton case. In the case of the KvBLL caloron with 

widely separated constituents one has to take the temperature scale, r = 47rTe~^^. Thus, the 2-loop recipe is to 
replace the factor {Sir^/g^)'^ [ke"^ / ^■nTf'^/^ in eqs.(80,83) by the r.h.s. of eq.(87). 

In contrast to the zero-temperature instanton, in the KvBLL caloron case this replacement is not the only effect of 
two loops. In particular, the potential energy P(v) is modified in 2 loops [33]. Nevertheless, the above modification is 
a very important effect of two loops, which needs to be taken into account if one wants to make a realistic estimate of 
the density of calorons with non-trivial holonomy at a given temperature. We remark that the additional large factor 
47re~'''^ « 7.05551 makes the running coupling numerically small even at T ~ A (l//3(r) ~ 0.07), which may justify 
the use of semiclassical methods at temperatures around the phase transition. This numerically large scale is not 
accidental but originates from the fact that it is the Matsubara frequency 2'kT rather that T itself which serves as a 
scale in all temperature-related problems. The additional order-of-unity factor 2e^''^ is specific for the Pauli-Villars 
regularization scheme used. 



VII. CALORON DENSITY AND INSTABILITY OF THE TRIVIAL HOLONOMY 



Since the caloron field has a constant component at spatial infinity, it is strongly suppressed by the potential 
energy P(v), unless v = 0, 2'kT corresponding to trivial holonomy. Nevertheless, one may ask if the free energy of an 
ensemble of calorons can override this perturbative potential. We make below a crude estimate of the free energy of 

non-interacting KvBLL calorons. We shall consider only the case of small v < ttT ^1 — --i=^ . If v exceeds this value 

the integral over dyon separations in eq.(80) diverges, meaning that calorons with holonomy far from trivial "ionize" 
into separate dyons. We shall not consider this case here but restrict ourselves to small values of v where the integral 
over the separation between dyon constituents converges, such that one can assume that KvBLL calorons are in the 
"atomic" phase. Integrating over the separation ri2 in eq.(80) gives the "fugacity" of calorons: 



c 

f{T/K) 



T^fiT/A) 



Sn^CjS'^ exp 



^_34\ 34 ln/3 



22 47rT 

P = -;r In n ) 

3 AeTJ^ 



nOO 

I{v) = / dRR^l + 2Tru{l-u)R]{2m^R+l)i''-^{2TT{l-iy)R+l)i^^-'''>-^exp 
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-2k R [--Av+Av^ 



(90) 



(91) 
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FIG. 5: Free energy of the caloron gas in units of T^V at T = 1.3A (dotted), T = 1.125A (solid) and T = 1.05A (dashed) as 
function of the asymptotic value of Ai in units of 2-kT. 



where we have introduced the dimensionless separation, R = ri2T, and the dimensionless v = 277^- One should be 
cautioned that eq.(80) has been derived for i? ^ 1, therefore the caloron fugacity is evaluated accurately if the integral 
(91) is saturated in the large- i? region. 

Assuming the Yang-Mills partition function is governed by a non-interacting gas of calorons and iV_ an- 
ticalorons, one writes their grand canonical partition function as 



Zeal = exp 



3 ^ ' 



E WJW-Ar''^) =eM~VT^H^,T)\, (92) 



where !F{v, T) is the free energy of the caloron gas, including the perturbative potential energy: 

47r2 

Hv, T) = - vf 2f{T/K) (93) 

We plot the free energy as function of v in Fig. 5 at several temperatures. The function f{T/A) rapidly drops with 
increasing temperature. Therefore, at high temperatures the perturbative potential energy prevails, and the minimal 
free energy corresponds to trivial holonomy. However, at T ~ A the caloron fugacity becomes sizable, and an opposite 
trend is observed. In this model, Tc = 1.125 A is the critical temperature where the trivial holonomy becomes an 
unstable point, and the system rolls towards large values of v where the present approach fails since at large v calorons 
anyhow have to "ionize" into separate dyons. 

Although several simplifying assumptions have been made in this derivation, it may indicate the instability of the 
trivial holonomy at temperatures below some critical one related to A. 
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APPENDIX A: ADHM CONSTRUCTION FOR THE BPS DYON AND THE KVBLL CALORON 

1. General ADHM construction 

The basic object in the ADHM construction [28] is the A; x (fc + 1) quaternionic-valued matrix A which is taken to 
be linear in the space-time variable x: 

A{x)^A + Bx, x = x^a^, ct^ = (if, I2) . (Al) 

The ADHM gauge potential is given by 

A^{x)^v\x)d^v{x), (A2) 
where v{x) is a (fc -I- 1) dimensional quaternionic vector, the normalized solution to 

A^(x)i;(x)=0, (A3) 
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and k is the topological charge of the gauge field. An important property of the ADHM construction is that the 
operator A\x)A{x) is a real- valued matrix: 

/ = (A(x)^A(x))-i G 11'=^'= . (A4) 

In what follows we shall use the equation 

A/At ^ 1 _ ^^t _ (A5) 

It becomes obvious when one notes that both sides arc projectors onto the space orthogonal to the vector v, which 
follows from v^v — 1, A^v — 0. 

In the case of finite temperatures, because of the infinite number of copies of space in the compact direction, the 
topological charge k = oo, and it is convenient to make a discrete Fourier transformation with respect to the infinite 
range indices. The Fourier transformed are 2 x2 matrix- valued functions t;(a;^, 0) of a new variable .2 € [—1/2,1/2] 
and A becomes a differential operator in z. 



2. ADHM construction for the BPS dyon 

As stated above, at nonzero temperatures the essence of the ADHM construction is the introduction of 2 x 2 
matrix- valued functions v{Xf^, z). The scalar product is defined as 

.1/2 

{vi\v2) = v^{x^,z)v2{x^,z)dz . (A6) 
J-1/2 

For the BPS dyon solution v has been found by Nahm [25] : 



■u(a;^,z) = . /— - — -exp{izvx^) (A7) 
y smh{vr) 

where aj^ = (I2, —ir), x^ = X/^iaj^ and r = \x\). The matrix-valued function v is the solution of the equation 

At {x)v{x, z) = 0, At (x) = id^ + vx^ (A8) 

normalized to unity, 

{v\v) = 1. (A9) 

The gauge field is expressed through v as 

= {v\d^v). (AlO) 

We use anti-hermitian fields such that the covariant derivative is = 3^ -|-^^. Comparing eq.(A8) with the general 
eq.(Al) we conclude that in this case 

A = -id,, B = v. (All) 

Eq.(A7) corresponds to the 'hedgehog' gauge. However we find it more convenient to work in the 'stringy' gauge 
where has a pure gauge string-like singularity. One proceeds from the 'hedgehog' gauge to the 'stringy' gauge 
using the singular gauge transformation (see e.g. [9]) 

v^v^ = vSl, A^-y A'^ = S-A^sl + S-d^^sl (A12) 

with 

S. = e-'^'e'^-^e-^^-' (A13) 
having the property that it "gauge-combs" A^ at spatial infinity to a fixed (third) direction: 

S-n^T^Sl = . (A14) 

In the 'stringy' gauge 

' VT o 

■ , , , exp[zv(ia;4 + rr"^)] . (A15) 
smh(vr) 

One can check that — (u^lS^w*) gives the M dyon field in the 'stringy' gauge as in eq.(13). We note that in the 
'stringy' gauge has a remarkable property 

v'{xi+n,x) =e''^^''v'{xi,x) . (A16) 
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3. ADHM construction for the KvBLL caloron 



Unfortunately, the original paper [4] does not present an explicit expression for v, the main ingredient of the 
ADHM construction. Wc could have used rcf. [5] but it seems that ref. [4] is more informative in some other respects. 
Therefore, we have to calculate v ourselves. 

From the point of view of the original ADHM construction t; is a quaternionic vector of infinite length since finite- 
temperature field configuration can be viewed as an infinite set of equal strips, the total topological charge in 
being infinite. The bracket is formally defined as a contraction along this infinite-dimension side; 

{v\v) = v^'v . (A17) 

The gauge potential results from 

A^{x) = v^x)d^v{x), D^ = di,+A^. (A18) 
The vector v{x) is the normalized solution of the equation 

A\x)v{x)=0, A{x)=(^^^_^y (A19) 

where B is a square quaternionic matrix, A is an (infinite) quaternionic vector, x = x^a^, = {it, 12). Introducing 
the notations 

v{x) = ^-Hx)(^-l^y u{x) = {B^-x^)-'X^ , (A20) 

eq.(A19) becomes 

{B^i -x^)u{x) = XK (A21) 

The inverse of the matrix (B^ — .t^) exists almost in all points. The points where it does not exists, are monopole 
positions. We are interested in those singular points that lie in the interval < 2:4 < 1 (we have rescaled the units to 
set temperature T = 1). The unknown function ^(x) is determined from the normalization of v: 

^x) = l+u^u. (A22) 

The formalism of infinite-dimensional matrices is not convenient. Following Nahm [25] we pass to the Fourier 
transforms in the discrete but infinite-range indices and get instead a continuous variable 2; e [— |, 5]. In the notations 
of ref. [4]: 



(fit _ a;t)„„ = - / e-2--"i^t (^) ^2.^zm^ (^) ^ + ^t(^) + 2mx^ = -d, + 2mr^{z) , 



-1/2- 

where 

A{z) = 2ni[^ + ■ aQ^{z)]. (A23) 

Here the function = 2a; when z G [—ui,ui] and —2u) otherwise; r^{z) = rfj,{z)aj^, x^ = x^aj^. As it can be seen 

from eq.(A23), the quaternion ^ simply represents the center of mass position of the whole system and can be set to 
zero, ^ = 0. We define r^(z) = when z £ [w, 1 — a;], and r^{z) — otherwise, where 

s = x — 2a)r72, f—x + 2ojfi2, S4 — = X4. (A24) 

Here u = \ — u), and and have the meaning of the vectors from the dyon centers to the 'observation' point, 
ri2 = r — s has the meaning of dyon the separation. We choose dyons to be separated in the 3d direction: ri2 = ri2e3. 
As for A, 

Xl = pC/te-2Tmd;.r^ (A25) 
where f/ is a unitary matrix, and p > 0. We have an additional constraint [4] 

^ . a = 7rp2 (i/^uJ • aU) ju) . (A26) 
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Here 2TTiu} ■ f is the value of at spatial infinity, lo = \lo\. It can be seen that np^ = ri2. We choose to rotate the 
direction in color space instead of rotating monopole positions, so we do not loose the generality of the solution. 
We connect the vector w and U by 

. nj = 0JT3 . (A27) 
Writing down the m}^ component of the (infinite) quaternionic vector as a Fourier transform 

Umix)= / u{x,z)e-^'''"''dz, (A28) 

J-l/2 

eq.(A19) we have to solve can be rewritten as 

{d^ - 2TTir^{z))u{x, z) = 2TTipU^ J2 = 2nipU^P+6{z - u) + P-5{z + u)), (A29) 

n 

where 

P± = ^{l±oj-T/uj) . (A30) 

Eq.(A29) is piece- wise homogeneous, therefore we present its solution in the form 

, , J exp(27rzs^z)i3i, —u<z<u f\o-]\ 

exp(27rirt(z-l/2))B2, iO<z<l-u; ^^'^^> 

and match the values and the derivatives of u at the endpoints of the pieces. 



where 



/i = 2wipU^P+ , /2 = 27ripU^P- , u>=--ij. 
Note that i?i,2 are matrices that generally do not commute: 

B2 = (e-'^^'^^'^e-^'''''^^ - e^'^'^^'^e'^''"'''^^ (^e"^''"*'^/! + e^""'^^ = 622&2ie"^'''''*'^^^f/VV', (A33) 

-1 



where 



622 = — cos(7ra;4)(chi shi f + chi shi s) + isin(7ra;4)(chi chi + f sshi shi) 

[_ 22 22 22 22 

612 = — cos(7ra;4)(chi shi f + chi shi s) + i sin(7ra;4)(chi chi + If shi shi) 
621 = 2nip{chi - srs shi), 611 = 2mp{chi + f T3 shi )e'''''*'^^, 

ib = — €08(2^X4) +chi chi H shi shi . (A34) 

Hat over the variable (notation found also in [4]) means contraction of the corresponding normalized vector with Pauli 
matrices, e.g. Cj = Q ■ t juj . We denote for brevity 

sh = sinh(47rsa;), ch = cosh(47rsa;), sh = sinh(47rra;), ch = cosh(47rra;), (A35) 

and the hyperbolic functions with subscript "i" arc the corresponding functions of half the same arguments. Com- 
bining eqs.(A31,A33,A34) back into eq.(A20) one gets the two-dimensional quaternionic vector v(x,z) which is the 
base for the construction of the Green's function, see Appendix B. Note that we have made a Fourier transform of u 
(A28) and got a continuous index z, so that scalar products of infinite-dimensional vectors become z integrations, see 
eq.(A38). 
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We note that U is actually a gauge transformation of v. Therefore, the gauge potential is obtained by a global 
gauge transformation of A'j^^^. We conclude that the determinant does not depend on the relative 'color orientation' 
of the Polyakov line or holonomy, and of the vector ri2 connecting monopole centers. Thus, we setU = 1 and lo = 0)63. 

We notice further that v{x, z) built above gives A^^ that is not periodic in time direction and zero A4 at spatial 
infinity. It is a peculiar feature of the 'algebraic' gauge used in [4]. It is more convenient to use the gauge in which 
the fields are periodic. To that end we make a non-periodic gauge transformation g = e^'^*^*'^'^^ and obtain 

v{x,zr^ = (^(^f^)) , w = u9, (A36) 

meaning 

w{x,z) =u{x,z) e^'^*^*'^-^^ , (A37) 
In terms of the Fourier-transformed v the bracket takes the form 

{v\v)=vlvi+ / vlv2dz, (A38) 
J-1/2 

where Ui is an upper element and V2 is a lower one. 

Now let us determine We use the following identities: 

b\^bi2 = ^2^22 = 4i^2i = 47rV^ (chi - y shi ) , b\ibn = ^t^'^P^ (du + ^ Su) • (A39) 

Note that the right-hand sides of eq.(A39) are proportional to the unity 2x2 matrix. Now we can easily calculate 
the normalization: 

{v\v) = ^-^b\^bnb\^bi2 f d2e-4-?---+^-25t^/,2i42&22 / e'^^^^'^^ 

tp — -ip 



2 

Tip 



chi shi chi shi \ r-n 

2- + + —shi shi 

s r sr ^ 



We used the identity f—s = ru = 1^1263. Thus for $ we get 



, / chi sh 1 chi shi \ ^2 

$=^, ^ = ^ + ri2\— 2. + ^ 1 +^shishi. (A40) 

^ I s r J sr ^ ^ ^ 

We have checked the of the KvBLL caloron (7) by calculating {v^'^'^ld^v^^'^). Note that v^'^'^ has the following 
periodicity property (only for integer n): 

v^P^'iz, X4 + n)= e2'^*"^t;P«'(^, X4) • (A41) 

APPENDIX B: SPIN-0 ISOSPIN-1 PROPAGATOR 

1. General construction of the Green function 

Once the self-dual field is found in terms of the ADHM construction, such that the gauge field is written as 
A^ — {v\dfj,v) where the scalar product is defined in cq.(A38), it is possible to construct explicitly the Green function 
of spin-0 isospin-1 field in the background of the self-dual field [4, 5, 25]. The solution of the equation 

- {DlY" {x)G'^'{x, y) = 5^' 6^'\x - y) (Bl) 

is given by 

^ab, ^ _ \T^T-{v{x)\v{y))T^{v{y)\v{x)) 
^""'^^ " A^^{x-yf 

+ / dzidz2dz3dziM{zuZ2,zs,Zi)-Ti{V\x,z^)V{x,Z2)T'')Tv{V\y,Zi)V{y,zs)T^) , (B2) 
47r j_i/2 ^ 



25 



where V(.T, z) = B^v{x, z). 

We denote the first term by G\ and the second term (the M-part) by G^- The only new object is the function 
M{zi,Z2, Z3, Z4) which we determine bellow. As we shall see, we do not need M with arbitrary arguments, but only 
at Z3 = Z4. For coincident arguments we obtain 

M{zi,Z2, z, z) = S{zi - Z2)M{zi, z), (B3) 

see below. 

The propagator (B2) is written for the IR'* space and does not obey the periodicity condition. The periodic 
propagator, however, can be easily obtained from it by a standard procedure: 

+•00 

G{x,y)= ^ G{x4,x;y4 + n,y). (B4) 

n— — 00 

In what follows it will be convenient to split it into three parts: 

e^^GiUo, e'-^^Gi, g^^J2^2. (B5) 

The vacuum current (20) will be also split into three parts, in accordance to which part of the periodic propagator 
(B5) is used to calculate it: 

j^^ = Ji.+j;+j:r- (B6) 



2. Propagator in the BPS dyon background 

In Appendix A. 2 we have found the needed periodic quaternion v{x, z) for the single BPS monopole (see eq.(A15)). 
The 4-argument function M for the BPS monopole was computed in ref. [25]. The result with the two last arguments 
taken equal is 

M(z3, Z4, z, z) = (5(z3 - Zi)M{z-i, z), 

M{z, z') = —^ {2\z -z'\-l + Azz') . (B7) 

Eqs.(A15,B7) completely determine the periodic propagator defined in eqs.(B2,B4) in the BPS dyon background. The 
use of this propagator is demonstrated in Appendix C. 



3. Propagator in tlie KvBLL caloron background 



In Appendix A. 3 we have found the needed quaternion v{x, z) for the KvBLL caloron. In this Appendix we derive 
the M-function for the KvBLL caloron. The propagator (B2) will be then completely determined in the caloron 
background. 

In the notations of [25] M is an infinite-dimensional rank-4 tensor, with indices nmning from 1 to k, the topological 
charge in IR^. As in the case of v, it is convenient to make the Fourier transformation with respect to the indices: 

/•1/2 

M„= / M{zuZ2,Z3,Zi)e''^'^-^''^'-^'^'*Uzidz2dz3dzi. (B8) 
The tensor Mpg^n is defined by the equation [26] 

^Tr[{AU)u{B^B)mj + {B^B)u{AU)mj - 2{A^ B)u{BU)mj]Mrsij = 6rl5sm , (B9) 

All indices here run from 1 to A; as rectangular k x (fc + 1) matrices A and B are contracted along the longer side. 
Here A and B are: 



A{x) = A + Bx, ^ = A(0), B = 




(BIO) 
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Eq.(B9) can be rewritten as 

iTr[(At A(0)),,(5„j + {A^AiO))mj5ii - 2BlB^j]Mrsij = 5ri5sm , (Bll) 

where A{0) ^ \ + B, B and A are found in eq.(A23) and eq.(A25), respectively. In our case k is infinite and 
we rewrite eq.(Bll) in the Fourier basis: 



1 

2^ 



At A(0, zs) + At A(0, + 2 ( ^ + rt (^3)) - r{z^) 



M{zi,Z2,Z3,Z4) = 6{zi - Z3)S{Z2 - Z4) , 



where r{z) = r^ai when z € [oj,! — oj], and r{z) = SiCi otherwise; ctj = ZTj. Zero components of r^, are absent 
because a;^ = 0. We use 

A^ A{Q,z) =-^+r^{z) + ^ {6{z-oj)+6{z+u)) . (B12) 

Here the first two terms come from the Fourier transformation of B^ B (A23) and the last one comes from the Fourier 
transformation of At A. We obtain the explicit equation for M: 



+ \ f{Z3) - f{Z4)f ) M{zi, Z2, Z3, Zi) + 



47r2 

+ "2" + 5(z3+uj) + 6{zi-u}) + (5(z4+w)) M(zi, 2:2,23, Z4) = 5(zi - Z3)5{z2 - Z4) . (B13) 

In the case Z3 = Z4, which is the only one we need as we shall see in a moment, we look for the solution in the form 

M{zi, Z2, z, z) = 5{zi - Z2)M{zi,z) . (B14) 
The equation for the two-argument function simplifies to 

where ri2 = Trp^. We see that the solution has to be piece- wise linear in its arguments. The solution is symmetric in 
its two arguments and for z < z' is given by 

32r^2'^^ui{z-ui){l-z'-uj)-8ri2-7T^{ui^+z{z'-l))+Tr ^/ c [/ , 1 _ , ,1 

2d(8ri2™ai+l) ' Z,Z <^ [L0,1 UJ\ 

M{z,z')={ ,e[_a;,a;],z'e[a;,l-a;] . (B16) 



2 \ Sri27rix>uj+1 ri 

4ri2 7r [z — z' — 2z'z—2{Lj — l)uj — 8ri2'^iz^—^){z+uj)u;]+'K 
2ri2(8ri27ra;w+l) 

Outside this range M is defined by periodicity: M{z + n, z' + m) = M(z, z'), where n, m arc integers. 

Now let us demonstrate that actually only the two-argument function M(z, z') is needed to construct the propagator 
satisfying the periodicity. It turns out that making the Green function periodic simplifies G2 (section III). One has 
from the definitions (B2)-(B5): 



^" = 1]^/ dzi...dz4M{zi...z^T^{y\x,z^)V{x,Z2)T")T:x(y\y^,Zi)V{y^,Z3y) , (B17) 
"-r" f^Ti" 7-1/2 



r-1/2 

( 

1/2 

where = 2/4 + n, = y. Using eq.(A41) we put 

V(y", z) = e2'^'''"^V(y, z) (B18) 

Further on, we note that for jryj < 1 one has 

Y,^{V\y^,Z4)V{y^,Z3)T^)=Tv{V\y,Z4)V{y,Z3)T^) 5(^3 - ^4) n • 
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Now we can see that making the Green's function periodic results in the substitution 

M{zi,Z2,Z3,Z4) -r^M{zi,Z2,Z3,Zs)S{zs - Z4) = ^M{zi, Z3)6{zi - Z2)S{Z3 - Z4) ■ 

\v\ \v\ 

It follows from eq.(A16) and eq.(A41) that for the monopole one has to take r] = v/(27r) < 1 and for the KvBLL 
caloron jj = 1. In both cases the M-part of the periodic propagator is given by 

^" = f4^ dzdz'M{z,z')TT{V\x,z)V{x,z)T^) Tr {V\y, z')Viy, z')t') , (B19) 

OT^ Wl J-1/2 

where the two-argument M functions are found in eq.(B7) and eq.(B16), respectively. 

APPENDIX C: VACUUM CURRENT IN THE BPS MONOPOLE BACKGROUND 

We compute the vacuum current (20) in the BPS monopole background in this Appendix. We assume < v < 27r 
and work in the stringy gauge (13) dropping the index s in given by eq.(A15). 

1. Singular part of the monopole current 

This part of the current corresponds to the second term in eq.(B5). At x ^ y this part of the propagator is 
singular. The regularization is presented in Appendix E. Eqs.(E2,All) state: 

J««^ = irtr(r'^j^), j^ = ^(v\fa^A^f\v)-h.c., AHx) = id, + vx^ . (CI) 
The function f{z, z',x) for the BPS monopole is known [25]: 

^ivX4(z—z ) , y^ . 

f{x; z, z') = (^sinh vs|2: — z'\+ coth — sinhvs2: sinhvs2;' — tanh — cosh vs^; coshvs2:'j . (C2) 

Here we denoted by s the distance to the M-monopole center. It is helpful to calculate the action of the Green function 
on V. Since monopole is a static configuration, we can take X4 = 0, moreover / is a scalar function and we can move 
Sl_ matrix to the left: 



, „ ^, y: cosh(si;2;) tanh(st;/2) — 22;sinh(sti2;) , smh(svz)coth(sv/2) — 2zcosh(svz) 

v) = S-f\v)\^^^^ = , / . . 12 + . / . . 

sv smn.\sv) A^ svsaA\{sv) 



1X4=0 

We use the following identities 



T3 • 



S-firTSl, = T3 , 

S-n$TSl_ = — cos(^)ti — sm{^)T2 , 

S-n^rSl = sm{(f))Ti - cos{(f))T2 (03) 



and arrive, after simple algebra, to 



{ii^jr-je,i4>) = T7r^i^\{^^ -T3,cos0 Ti+sin0 r2, -sint/) n+cosfj) T2}{d:^ - yt3.s)\v) + h.c. . 

iZTT 

Finally we obtain the singular part of the vacuum current: 
Jl = 0, 

iv (s^v^csch^(sv) + svcoth(sv) - 2) _ ^ ,^ ^ . / 

Jl = ^ OA 2 2 ■ ^.1 ^ ^ Ti cos +^2 sinO^ , 

^ 2iiT^s-^ smh(sv) 

Je = ^ OA 2 2 ■ u ^ ^ Tisin ^ -T2COS , 

247r^s^ smh(sv) 

i (1 - s^v^ coth(sv)csch^(sv)) 
= 2A^^ ' 



where (Tg)" = 
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2. Regular part of the monopole current 

We are going to calculate the part of the current that corresponds to 

igT\x,y) = ^ r, 2(\ ^2 Tr [T''{v{xMyn))T''{v{ynMx))] , yn = y, ynA = yi + n, 

yn) 



(C5) 



namely 



At first consider J^^. We have to compute with equal arguments. Substituting (A15) into (C5) and calculating 
the trace one has 



Tr 



(C6) 



where 



„2 2 inv(z-z') 

Tr = 2 _ [cosh(2sv(^ - z')){S''^ - S'^^S^^) + cos]i{2sv{z+z'))S^Sl + smh{2sy{z-z'))Tf] 



sinh (sv) 



with T"^ = is°'^°. To compute the sum in this expression we use the summation formula (note that v < 27r) 



e«n _ ^ _ jzl J_ 



-2tt < Z <2'K . 



(C7) 



It remains now to calculate integrals over z and z'. The result is 

,2 



3coth(sv) — sv(3csch^(sv) + 2) _^ (svcoth(sv) — 1)^ 



Sirs 



sv cs(ii^(.s\-) — (•c)ili(.sv) 1 — (•S(ii"(,sv) 1 



Stts 



167r2s2 



12 



5"^ . 



Now we turn to the J^^ part of the current where wc have to sum over n a derivative of the propagator. First of 
all we consider derivatives of the trace in (C5). One finds for x = y: 



2s^y^i 



ivn{z—z') 



{d^ - d^)Tr= . (cosh(25v-g) + cosh(2sv^'))(I'i sin - T2 cos0)e 

SlUJll 5V ) 

(55 - 5^)Tr= fsin6> J \^ cosh(sv(z + z')) cosh(sv(z - z'))(Ti cos<p + T2 sine 
1' V \^ smh(svj^ 



cosh(2sv(^ - z')){l - cose)Ts ] e 



sinh(sv)' 

.23- 

{dl - 9|)Tr = - . ' {z - z') sinh(2sv(0 - z'))ne'^''^^-^'^ , 



Jvn{z-z') 



sinh (sv) 



(5,--S,^)Tr = 0. 



(C8) 



Here only terms even in ^; — were left. The last two equations are especially clear as we can drop out the matrices 
S in eq.(A15). 

A derivative of the denominator of (C5) is equal to zero for a; = y except for the derivative with respect to X4, but 
in this case wc have the expression of the form (C6) with instead of in the denominator. Now we can sum 
over n. We use the summation formula 



Ee"^" _ z'' z\z\ z 



67r2 27r 3' 



-2w < z <2tt . 



(C9) 
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Next one has to integrate over z, z' . Combining all pieces we obtain: 

Jl = 0, 

_ . cos((/))ri +sin(0)r2 
^ ' 48sinh3(sv)7r2s2 
^ .sin{cj))Ti - cos{<j))T2 
48 sinh'*(sv)7r2s2 

where we denote 

pi = (s^v^ + Otts^v^ + 3v + 3s(v + tt) sinh(2sv)v - (s^v^ + 3v + Qn) cosh(2sv) + Gn), 
p2 = 87r^.s^(-l + svcoth(sv)) - 127rscoth(sv)(-l + svcoth(sv))^ 

+ (-3(1 + 4s2v^) + sv(4(3 + s^v^) coth(sv) + 3sv(-4 + svcoth(sv))csch2(sv))). 

Wc have used spherical coordinates. For example, a projection of J onto the direction ne = 
(cos 9 cos 0, cos sin (/>, — sin 6*) is denoted by Je . 

3. M-part of the monopole current J™ 

Combining together eqs.(B7,B19) we have for the M-part of the periodic Green's function: 

V /"^/^ 

0"^ = -— / dzdz'{2\z-z'\-l + 4zz') Ti{v\x,z)v{x,z)t'') Ti{v^{y,z')v{y,z')T') (Cll) 

Note that we can drop out S- in (A15). In the stringy gauge one has 

v^{xij,,z)v{x^,z) = ■ w X exp[-2vsT^^] . (C12) 

SlUill V5 1 

It means that ^™ has only the '33' component. Taking the trace we get: 

Tr [v^{x^,z)v{x^,z)t^] = -2 . sinh(2vs^;) . (C13) 

Therefore the only nonzero component of G™ is 



.33.. _ rdzdz'{2\z-z'\-l + 4zz') ^inh(2vr,.0sinh(2vr z) 

47r J -1/2 smh(vrx)smh(vry) 



{Vx = \x\, ry = \y\). Performing the integrations we get 



33 1 / 1 , vr^cothvr^ -vr^cothvr^ 

ix,y) = ^— + ) ■ 

Note that (C14) is symmetric in its arguments. For that reason the contribution to the current coming from (ordinary) 
derivatives of 0™ is zero, 

dxg"'{x,y)-dyg'^{x,y)\^^^ = 0, 
and only the anticommutator x), A^} remains. Taking the limit x — > y wc get for ^™ 



w X 2 VS 

coth(vs) h ^ 2 



VS sinh (vs) 



(C15) 
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and for the contribution to the current J™ = A^Q"^ + Q'^An we obtain in spherical coordinates 

J?" = 0, 

_ i(sin((;i)T2+cos(0)Ti)(svcoth(sv)+s^v^csch^(sv)-2) 
■^^ ~ IBtt sinhr.svVs^ ' 

7m . 
Jo ■ 



IGtt sinh(5v)5^ 

i(siii((;i)Ti-cos((;i)T2)(sv coth(sv)+s^v^csch^(sv)-2) 
167r sinh(sv)s'^ ' 



Jf = 0. 



(C16) 



Adding up (C4,C10) and (C16) we obtain the full vacuum current in the BPS background, see eq.(23) of the main 
text. 



APPENDIX D: VACUUM CURRENT IN THE KVBLL CALORON BACKGROUND 



There are no principal problems to make the calculation of the caloron Green's function and the ensuing vacuum 
currents exactly. One can consider this Appendix as an instruction how to perform the exact calculation. In fact, 
we have done it but unfortunately the exact result for the current is about 200 pages long and thus too large to be 
printed. However, in certain limits the expressions drastically simplify. In particular, assuming the case when the 
dyons inside the caloron arc widely separated such that their cores do not overlap, it is relatively easy to find the 
KvBLL caloron current with the exponential precision (i.e. dropping out term of the order e"*^^, e""^). This will be 
sufficient to find the determinant of the KvBLL caloron for large r\2 up to some constant. 

With the exponential precision, the only nonzero components of the KvBLL caloron's gauge potential in fundamental 
representation are (see section II) 



IT'i ( 1 



ir3 fl 
' 2 [r 



{ri2 - r + s)(ri2 + r - s) 



2 \ " ' r s J ' ^ 2 \r ' s J V (ri2 + r + s)(r + s — ri2) 
We are using the coordinates X4,r,s,(p, where r, s are defined in (A24) and is defined by 



X = Xg{cos(p eb + sintp ei)+ 
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2n, 2n2U.je3, X,- 

One can easily check the consistency of this definition, i.e. that 

s = r = |f|, where s = x — 2u!ri2, f=x + 2uiri2, n3 = ?"i2e3 . 

Since Ar — and As = we have to calculate only the J4 and components. 

We shall use the ADHM construction. The main steps of the calculation are similar to that for the monopole. 
Dropping out exponentially small terms in eq.(A36) one has in the periodic gauge 



_ \^{ri2+r-s){ri2+s-r){r+s-ri2){ri2+r+s) 



(Dl) 



• (D2) 



v{x, z) 



r + s + ri2 V w{x,z) 



(D3) 



w{x, z) 



-iirpe 



27r(is^ z—LUs) 



ri2+i — s 



T3 



2tt {ir^ z — ir^ /2—rLj) 



(ri sm.^p + T2 cos<^) 



(il2 sin7ra;4 + ra cos7ra;4) + 



r+s-ri2 



2a; „ 



(sin(7rX4 — ^p)ti — cosinxi — ifi)T2) 



-U<Z<ljJ 



U)<Z<1—U) 



(D4) 



r+s—rri 

where = s^uL s = S4 = 0:4. We shall use the following formulas to pass to the cylindrical coordinates x^^Xg, 



dr 
dxs 



+ — 



2ri2r 



dr Xg 
dx„ r 



dr 
dip 



= 



ds 

dX3 



r\2 + — 



2ri2S 



ds Xg 

dx„ s ' 



5s 
d(p 



= 



(D5) 
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1. Singular part of the caloron current 

Let us calculate the singular part of the vacuum current with exponential precision. It is related to the zero 

Matsubara frequency. Similar to the monopolc case, wc could use cq.(E2), whore the Green's fmiction (A4) for the 
case of KvBLL caloron was found in [4]. However it is more convenient to use eq.(ElO)) because then we have only 
to take derivatives of the simple expression (D3) and no integrations arise. Eq.(E2) would have been more suitable 
for the exact calculation. 

It is straightforward to calculate the quantity T from eq.(E9). It is sufficient to calculate the second time derivative: 



r ~< d^divlv > -Al . 



(D6) 



Bearing in mind that F is a vector under gauge transformations, we can perform calculations in any gauge. Up to the 
exponentially small terms we have 



po6 



(r + s) ((r - sf + ri2 (r + s)) 5"^ 
4^2 g2 ^^^2 + r + s) 



(D7) 



One can observe from eq.(Dl) that all terms with derivatives in the right-hand side of eq.(ElO) are zero. Writing the 
Laplace operator in the cylindrical coordinates we find 



1 



487r2 
1 

487r2 



{dJdi+di+—dgXgdg+^d'' ]v\v)-h.c 



247r2 



{Al+A^AiA^+&AiT) 



+ ^^{Al+A^A^Ai^-QA^T) . 



Taking the derivatives we obtain simple expressions: 



31 = 



34, 



ITS 

487r2 



1 



iT3Xgri2 



1 1 

r s J STr^rs{ri2 + r + s)^ 



(D8) 
(D9) 



2. Regular part of the caloron current 

Next we calculate the temperature-dependent part of the KvBLL caloron vacuum current. As in the monopole case 
(Appendix C.2) we divide the current into two parts, 

Jl = Jf + Jl', (DIO) 

where 

jr2 _ V- 1 ^^, Tv[T-v^x)v{yyv\y)v{x)] 



Tr[T''vT{x)v{y)T^vT{y)v{x)] 



and 2/4 = X4 + n. The quaternion function v{x,z) has been constructed in Appendix A. 3 (actually called v^^'^{x,z) 
there). It is important that v{x, z) has the remarkable periodicity property (A41). 

In evaluating the above currents the tactics is to factor the matrix part out of the integrals over z. We use the 
following notations for the integrals over z: 

e^^^*"^ cosh(47rsz)dz, 11= e^'''''^'+^/^'^ cosh{A'Krz)dz, 

-UJ J — U) 

r = r e^'^'"^ sinh(47rs^)d^, = C e^'^'^^^+^s) sm\i{4-Krz)dz . 
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We obtain the following relations for the matrix structures: 

i,''wHxi)w{x4. + n) = {IIP + + {11^ + IHPr), 
i>^w^{x4 + n)w{x4) = {1113 - I1f3s) + {11^ - Il^r), 
iP'^w''{xi)d4w{xi + n) = {I^Po + I!1P° + idsllfi + idJlP,) + ... 
i;^w\x4, + n)diw{xi) = {Ill3o - - idJUP + idsIlPs) + ... 

where means the same expression but with bar over each quantity and r instead of ,s. The notation means: 
derivative from the right minus derivative from the left. The definition and the evaluation of the matrix structures 
with the exponential precision is 

4° = blbl2rd4b22b2i ^ '-^{^ + 1)(S3 - l)e«-^'^e4--^, 
/Jo = b\,b\2dibi2bn ~ o{e^^''^e'^n, 

2 

Po = 4i425462262i ^ ^(i? + 1)(S3 - l)we8-«-e^'^™, 

= 6I16I2S612611 ^ ^(i? + l)(r3 + l)we*-«-e8-™, 
Pr = bl,bl2fb22b2i ^ ^(t? + 1)(S3 - l)we«"*'^e*"™, 

f3 = 6li4&i26ii l)(r3 + l)e4--e«^'^^, 

= blbl2b22b2i c^^-^{^+l){l- ss)e^^'"e'^^^, 



where 



^ ~ + I)g4.«a,g4.™^ ,? = = ^, rr3 = SS3 - 



4 sr 2sr 2ri2 

Substituting this into in the currents JJ^^'^^ we obtain certain sums, which are of the form 

^/^/^/(47r2n2), ^cos(27rna))3J^/(47r2n2), ^ sin(47ma;)/(7r2n3). 



All such sums can be calculated using the summation formulae 

g2^m(z-l/2) ^2 1 , , 1 1 

^ 4n-n^ =y-24"^^("^' -2^^^2' 



For example, 



g2,rm(^-l/2) ^ X ^ j 



/'Cl.' /"•(J-' 
/ C2(^ + z'- 1/2) cosh(47rs^) cosh(47rs^') dz' 
-OJ OJ 



and so on. With some help from Mathematica we come to the final result 



r — 
4 



1 



1 



1 



4^2 j,3 ^2 j,2 ^ 



+ 



8 Stt 



4 TT^ .S-^ 



2 12 2 

irrs IT s'^ 3r 3s 



(Dll) 



nr'' nrs ns'' r s 3 



16 16 \ 2 647rw3 

16 77 I + 

r s 



in 
2 ' 



1 1 
+ 



iT3Xgri2 



r s J 47r2rs(ri2 + r + s)^ 



(D12) 
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3. M-part of the caloron current J™ 

This part of the current is especially simple: with exponential precision it is zero. The main steps are the same 
as in the case of a single monopolc. The starting formula is our eq.(B19). First of all we note that only the lower 
components of v are left and only the a = 3 component is nonzero; 

Tr [V+{x,z)V{x,z)t''] = ^Tr [w+{x,z)w{x,z)t''] oc . 

Inspecting the definition of the M-part of the propagator (B19) we observe that 

G'^''''{x,y) cx (5^3,5''^ G'^''''{x,y) = G""°^(?y,x) . (D13) 

The second equation means that the terms with derivatives in the expression for the current (20) cancel each other. 
It follows from the first one that the product of G™ and oc e^"^ is equal to zero, too. Therefore we conclude that 

j;' ~ . (D14) 



APPENDIX E: REGULARIZATION OF THE CURRENT 



Here we consider in more detail J^, the contribution to the current from the singular (as x y) part of the 
propagator Q^{x,y) defined by eq.(B5). This part is obviously temperature- independent, so the zero-temperature 
results are applicable. We regularize the current by setting x — y = e and inserting a parallel transporter to support 
gauge invariance, see e.g. [20]: 

TS TSl _i_ TS2 

Jf = [A^{z - e/2)G\z - e/2, z + e/2) + G\z - e/2, z + e/2)A^{z + e/2)] Pexp A^dz^^ , 

= [id; - y)] Pexp A^dz^^ , (El) 

where x = z — e/2, y = z + e/2 and wc imply averaging over all directions of e in the 4d space. This regularization 
method was proved to be equivalent to the ^-function regularization approach [29]. 

For a background field written in terms of the ADHM construction, a useful expression for the vacuum current was 
derived in refs. [20, 29] . In the SU (2) case it acquires the form: 

= ^^"''tr {r"j,) , j, = Y^ Wf (a^AtB - B^A^l) fB^\ v) 

(see Appendix A for notations of the ADHM construction elements). 

We would like to derive another expression for this part of the current - in terms of derivatives. In some cases it is 
more useful. We start from writing our result: 

^'^ = [(^.^'(H) 1^) - h.c.] . (E2) 

Let us prove it. First of all we consider the action of one derivative 

D^{vix)\=dM - df^{v\v){v\ = - \v}{v\) 

= d^{v\AfA^ = -{v\d^AfA^ (E3) 

= -{v\Ba^,fA^ . (E4) 

At the end of the first line we have used eq.(A5). The first equation in the second line comes from differentiating the 
ADHM equation 

= d^{{v\A) = d^{v\A + {v\d^A. 
The last equation follows from the definition (BIO). Therefore we obtain 

D^{v{xMy)) = -{v{x)\Ba^f,Ai\v{y)) = -{vix)\Ba^Mx - y^B^Hy)) , 
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where in the last hne wc have used the ADHM equation (A3). We next consider two derivatives. It is important here 
that / is proportional to the unity 2x2 matrix. We have 

Dl{v{x)\v{y))=-D^^{{v{x)\Ba,f,{x-y)^BMv{y))) 

= {v{x)\Ba^UAlBa^Ux - y)^B^v{y)) - {v{x)\Ba^d^Ux - y)^B^\v{y)) - {v{x)\Ba^f,alB^\v{y)) 
= -A{v{x)\Bf,B^\v{y)) . 

We have used here 

(tM- = -/x9^KAt A)/, = -/,c7^(at et A + A^Ba^)U = -2/xBt A/, = /^a^Atfia^/, . 
We have also used that the derivative of the inverse operator is d{0~^) = —0~^{dO)0~^ , as well as the relations 

a^aj; = 4, C7^c cr^ = -2c^ (E5) 

where c is an arbitrary quaternion. 

Finally, let us consider three derivatives: 

D;Dl{v{x)\v{y)) = -ad; {{v{x)\Bf,B^v{y))) = A{vix)\BUa^AiBUB^\v{y)) + A{v{x)\BdJ,B''\v{y)) . 

Notice that the last term is hermitian at x = y. Thus we have proven that the current written in form of eq.(E2) is 
equivalent to that of eq.(E2): 

^ [{D^D^v\) \v) - h.c] =j^{v \Bf (a^AtS - B^Aal) fB^\v) . 

In fact it is more useful to rewrite everything in terms of ordinary rather than covariant derivatives: 

{D^,D^{v\) \v) = {d^,d\\v) + A^A^A^ - d.A^Af, - A^d^A^ - d^A^A^ + d^d.A, + 6A^r , (E6) 
where is in the fundamental representation and 

i [D^D, + D,D^{v\) \v) = S^^T . (E7) 

Wc have to prove that the left-hand-side of eq.(E7) is a Lorentz scalar as is the right-hand side. Note that eq.(E5) is 
proportional to x — y. The only way to obtain a nonzero result at x — y ^ is to differentiate this factor: 



^ {D^D, + D,D^{v\) \v) = -l{v\B{a^at + a.al)fB^\v) = -S^.{v\BfB^\v) . 
It follows from eq.(E8) that F is hermitian. We can write F as follows: 



(E8) 



r<5^. = {d^,d^v\v) + \{d^A, + d^A^) - \{A^A, + A^A^) . (E9) 
Finally, the regularized singular part of the current can be written as 

3^ = ^ {{d^.d^v\v) - h.c.) + ^ {A^A^A, + d^d^A, + 3A^r + ^TA^) . (ElO) 
Eqs.(E9,E10) are used for the calculation of the singular part of the vacuum current in Appendix D.l. 
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